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1. Introduction 


The field of linear accelerators has grown to be so voluminous that compressing 
it into a few elementary lectures, as is our intention, requires one to be selective in 
some arbitrary way. In this first section, we will list some of the broad topics which 
the subject of linear accelerators encompasses and then indicate which of these will be 
covered in some detail. _ 

b.l The Accelerating Structure and the Beam 

Linear accelerators consist of two basic elements: the accelerating structure and 
the particle beam. 

The accelerating structure depends on the type of linac. The main types of linacs 

are: 

1. DC linacs, like Van de Graafs, in which the structure consists of some kind 
of column of electrodes. These electrodes sustain a DC electric field which 
accelerates a continuous stream of particles. DC linacs are limited to a few tens 
of MeV. 

2. Induction linacs in which the accelerating electric fields are obtained, according 
to Faraday’s law, from changing magnetic fluxes. These changing magnetic 
fluxes are generated by large pulsed currents driven through linear arrays of 
magnetic toroids. The beam path, along which the electric field develops, can 
be considered as the single turn secondary of a transformer. Induction linacs 
are generally used in medium-energy high-current pulsed applications. 

3. RF linacs, the type on which we will concentrate here, can be categorized in 
a number of ways: low frequency (UHF), microwave frequency (L, S, C, or X- 
band), laser frequency; CW or pulsed; traveling-wave or standing-wave; room 
temperature or superconducting. In all these cases, the structure is a conduct¬ 
ing array of gaps, cavities or gratings along which rf waves with an electric field 
parallel to the beam can be supported and built up through some resonant pro¬ 
cess. RF linacs are used for a wide spectrum of applications from injectors into 
circular accelerators, to entire high-energy accelerators such as SLAC, medical 
accelerators, and many others. 
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While the accelerating structure can be considered the heart of each individual 
machine, it cannot work without its associated systems such as the power source, the 
vacuum, cooling, support and alignment, and instrumentation and control systems. It 
is important to note that a technological breakthrough in any one of these associated 
systems can have profound effects on the main design of the accelerator. 

The second basic element of the linac is, of course, the beam and the particles which 
compose it. The vast majority of linacs today are electron (or positron) machines; they 
number about 1400 accelerators of which a large percentage are commercial radiation 
therapy machines. The other linacs accelerate protons (// ± ) and in a few cases ions. 
There are about 50 proton or ion linacs in existence in the world. 

Very little will be said in these lectures about the sources of these particles, guns, 
ion sources, duoplasmatrons, polarized beams, strippers, positron radiators, etc. This 
is an encyclopedic subject by itself with many specialties and sub-specialties which 
cannot profitably be summarized here. The reader should not conclude that because 
the subject has been left out, it is not of crucial importance to the design and operation 
of a given linac. Not only does the source have an effect on how well the accelerator 
can perform its function, but in some cases it determines how a new concept can or 
cannot be approached. An example of this is the electron source for the linear collider 
for which a conventional electron gun cannot create a beam with a sufficiently small 
emittance and must be followed by a damping ring to “cool” it down. 

The fundamental problems in beam dynamics are: 

1. Longitudinal bunching and stability, 

2. Focusing and transverse stability, and 

3. Steering and transport to a target or interaction area. 

As long as the number of particles or current density is not too high, the analysis 
of many of these problems can be carried out by following single particle trajectories 
and then formulating the behavior of the beam in terms of its envelope. The important 
parameters are transverse emittance, bunch length and energy spectrum. 

Far more difficult are the beam dynamics associated with collective effects such as 
emittance growth, beam break-up, beam-wall and wakefield effects in general. These 
problems are treated in other chapters of this book. 
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1.2 Plan of the Lectures 

These lectures come in five sections. The first is this introduction. The second is 
a short chronology of what are viewed as important milestones in the field. 

The third covers proton linacs. It introduces elementary concepts such as transit 
time, shunt impedance, Q, but defers a more systematic discussion of periodic ac¬ 
celerating structures to the next section. Critical issues such ais phase stability and 
transverse forces are discussed. Finally, examples of existing facilities are given and 
modern developments are mentioned. 

The fourth section contains an elementary discussion of waveguide accelerating 
structures. It can be regarded as an introduction to some of the more advanced treat¬ 
ments of the subject, such as are contained in the book on Linear Accelerators, edited 
by P. Lapostolle and A. Septier in 1970, which, except for some recent developments, 
is undoubtedly the best single source of detailed information on the subject. The en¬ 
cyclopedic nature of this volume, to which we will henceforth refer to as L.A. , partially 
excuses the sparsity of other references given in these lectures. 

The final section is devoted to electron accelerators. Taking SLAC as an example, 
various topics are discussed such as structure design, choice of parameters, frequency 
optimization, beam current, emittance, bunch length and beam loading. Recent de¬ 
velopments and future challenges are mentioned briefly. 

Since proton and electron machines will be discussed separately, it may be useful, 
before starting, to take a broad look at the full spectrum of linear accelerator appli¬ 
cations shown in Fig. 1. In this figure, the kinetic energy Ek, which is the variable of 
greatest interest for most physics applications, is shown in the horizontal axis and the 
normalized velocity (in the form /? — v/c), which determines much of the linac design 
is plotted along the vertical axis. For particles of a given mass mo, these variables are 
related by the formulae: 


E k ~ 


mpc 2 


2 

moc 



mpv 2 (non — relativistically) 


( 1 . 1 ) 

( 1 . 2 ) 
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KINETIC ENERGY (E«) 


Fig. 1. The full spectrum of linear accelerator applications, displayed at appropriate 
energies and velocities. 














































The curves which result for the various particles are plotted on the figure and 
facilities of successively higher energies are indicated. It can be seen that electrons 
are almost always relativistic while protons (not to mention heavier ions) almost never 
are. It is for this reason that proton and electron accelerators tend to be described in 
different terms. 
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2. Some Historical Milestones 


The history of linear accelerators is briefly described by J. Blewett in the book L.A. 
A technical discussion of some early projects is contained in the “Linear Accelerator 
Issue” of the Review of Scientific Instruments (February 1955). 

Here we have chosen to enumerate only some of the important milestones. The 
choice has led them to be spread more or less uniformly-in time, emphasizing the 
remarkable continuity of progress over a period of more than half a century which 
promises to continue into the future. This uniform spread results perhaps in the 
slighting of some of the early work when parallel progress was being made at many 
laboratories. 


The history of linacs can be viewed as a sequence of attempts to fool charged 
particles so that they see cumulatively acting voltages across linear arrays of gaps. A 
chronology follows. References can be traced through L.A. 

1924 A theoretical paper by G. Ising, Stockholm, describes a method for accelerating 
positive ions (canal rays) by applying the electrical wavefront from a spark discharge 
to an array of drift tubes via transmission lines of successively greater lengths. 


1928 An experimental paper (including the theory of the betatron) by R. Wideroe, 
Switzerland, describes the successful acceleration of potassium ions to 50 kV. Figure 
2 indicates schematically the setup in which the ions pass successively through three 
drift tubes: the first and last are grounded, the center one is attached to a 1 MHz 
oscillator with a voltage of 25 kV. The distance d between gaps is adjusted so that 


d 

v 



Ao 

2 


where / is the frequency and Xq the free space wavelength at that frequency. The 
potassium ions travel from one gap to the next in one-half an RF period. Since higher 
frequency oscillators did not exist at the time, lighter particles traveling faster could 
not be accelerated. 
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Fig. 2. The Wideroe linac. 

1931-35 K. Kingdon (G.E.), L. Snoddy (Univ. of Virginia) et al., accelerate electrons 
from 28 keV to 2.5 MeV by applying progressive wavefronts to a drift tube array. 

1931-34 E. Lawrence, D. Sloan et al., (U.C. at Berkeley) build a Wideroe-type linac 
(i.e., d = /?Xo/2) with 30 drift tubes, oscillating to 42 kV, driven by 7 MHz oscillators. 
Mercury ions are accelerated to 1.26 MeV. Oscillators of high enough frequency for 
protons are still not available. Similar work continues at Cornell, in Japan and in 
England. 

1937-45 W. Hansen and the Varian brothers invent the klystron (at first a low 
power device) at Stanford. Subsequently, the high power magnetron (2 MW pulsed) is 
developed in Great Britain for radar purposes as part of the war effort. 

1945-47 L. Alvarez, W. Panofsky et al. (U. C. Berkeley) build a 32 MeV proton drift 
tube linac (Fig. 3) three-feet in diameter, forty feet long, powered by 200 MHz war 
surplus radar equipment. As indicated, the Alvarez structure differs from the Wideroe 
structure in that all tubes are contained in one large cylindrical tank and are powered 
at the same phase: the distance between drift tubes is arranged so that the particles 
are shielded from the fields when they are in the decelerating phase. Adequate beam 
acceptance required that the accelerating field not have much variation with radius, 
thus precluding operation at higher frequency. As will be explained later, longitudinal 
phase stability turned out to be satisfactory but transverse focusing was problematical. 
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Fig. 3. An Alvarez linac. 


32 MeV 
Proton 


1947-48 At Stanford, W. Hansen, E. Ginzton, W. Kennedy et al., build the Mark I 
disk-loaded linac yielding 4.5 MeV electrons in a nine-foot structure powered to 1 MW 
at 2856 MHz. It is the first of a series: Mark II (40 MeV); Mark III (1.2 GeV); and 
SLAC (30 GeV). Parallel efforts take place in Great Britain, France and the USSR, 
and at M.I.T. and Yale in the U.S.A. 


1952 J. Blewett at BNL shows that alternating-gradient focusing works with quadru- 
pole coils inserted in the drift tubes, solving the transverse focusing problem for 
protons. The Alvarez linac (with some modifications) from then on serves as the 
model for most subsequent proton and ion linacs (up to 200 MeV). 

1965-66 P. Wilson, A. Schwettman et al., at HEPL, Stanford, report the first suc¬ 
cessful operation of a superconducting linac producing 500 keV electrons with three 
lead-coated cells. 

1967-69 V. Sarantsev et al., at Dubna, USSR, build a linear induction accelerator, 
as do D. Keefe et al., at LBL, shortly thereafter, both groups with the intent of 
accelerating electron rings. 

1971 R. Koontz, G. Loew, and R. Miller at SLAC for the first time accelerate a single 
electron bunch through the 3-km linac and show experimentally that beam loading is 
energy independent. 

1972 D. Nagle, E. Knapp et al., at LASL, Los Alamos, successfully operate their 
800-MHz side-coupled cavity linac, LAMPF, and produce 800 MeV protons. 
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1973 P. Wilson, D. Farkas and H. Hogg at SLAC invent the rf energy compression 
scheme called SLED which in the next five subsequent years gets installed on the 3-km 
linac, boosting its energy up to 30 GeV. 

1980 R. Stokes et al., at LASL, successfully test RF quadrupoles (up to 2 MeV) 
following a 1970 suggestion of I. Kapchinskii and V. Teplyakov, ITEP, Moscow. 

1982 H. Grunder, F. Selph, et al., at LBL, use the HELAC-and Bevatron to accelerate 
{7 238 with charge state ~ 69. 
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3. Proton Linear Accelerators 


3.1 General Features 

As was mentioned previously, the early linear accelerators were restricted to pro¬ 
tons and heavier particles. Vacuum tubes with inferior high frequency characteristics 
and various discharge devices were the only RF power sources available. In a single cy¬ 
cle of a 1 MHz oscillator, a highly relativistic particle travels 300m - a length hard to 
conceive tor a piece of apparatus in the early days of accelerator physics. A frequency 
more typical of modern proton accelerators is 200 MHz. The corresponding free space 
wavelength Xq is 


Xo = j = 1.5 m 

If the input velocity of a proton is such that 


v 

- = 0.05 
c 

corresponding to a kinetic energy 


(3.1) 


Ej( — 1.17 MeV , 


then the distance it travels in one cycle is 

v 

d = - Xo = 7.5 cm . 
c 

If we use drift tubes to shield the proton from decelerating forces, their length must 
be at least half this distance. Precise parameters for the FNAL 200 MeV proton 
linac injector are given in Table 1. It can be seen that, even at the output end, 
the particles are not highly relativistic. This means that the much higher frequency 
periodic structures typical of electron accelerators are not yet appropriate (for reasons 
of inefficiency and radial effects to which we will return). 

3.2 Transit Time Factor 

Consider a proton passing through a pill-box cavity as shown in Fig. 4, in which the 
maximum potential difference is Vq. If the gap length is g and the instantaneous electric 
field E z in the gap is independent of the longitudinal coordinate z , the maximum 
electric field is given by 
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Table 1. FNAL 200-MeV Linear Accelerator Specifications 

Output energy 
Output momentum spread, 
total for 80% of beam 
Peak Beam current 
Emittance at 200 MeV 
(each transverse mode) 

Beam pulse lcr^fh 
Pube repetition late 
Cavity resonant frequency 
RF pube length, variable to 
RF duty factor, maximum 
Synchronous phase angle, 
from rf peak 


Cavity number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Proton energy in (MeV) 

0.75 

10.42 

37.54 

66.2 

92.6 

116.5 

139.0 

160.5 

181.0 

Proton energy out (MeV) 

10.42 

37.54 

66.18 

92.60 

116.5 

139.0 

160.5 

181.0 

200.3 

Cavity length (m) 

7.44 

19.02 

16.53 

16.68 

15.58 

15.54 

15.83 

15.88 

15.73 

Cavity diameter (cm) 

94 

90 

88 

88 

84 

84 

84 

84 

84 

Drift-tube diameter (cm) 

18 

16 

16 

16 

16 

16 

16 

16 


Bore-hole diameter (cm) 

2.0 

9 ^ 

3.0 

3.0 

3.0 

4.0 

4.0 

4.0 

4.0 

4.0 

Cell length L 










(first cell) (cm) 

6.04 

22.2 

41.1 

53.3 

61.8 

68.2 

73.3 

77.6 

81.3 

(last cell) (cm) 

21.8 

40.8 

53.8 

61.5 

67.9 

73.1 

77.4 

81.1 

84.3 

Gap length G 










(first cell) (cm) 

1.30 

4.4 

12.2 

19.5 

22.6 

27.1 

30.9 

34.3 

37.3 

(last cell) (cm) 

6.70 

12.7 

19.3 

25.1 

26.9 

30.8 

34.2 

37.1 

39.7 

G/L 










(first cell) 

0.21 

0.20 

0.30 

0.37 

0.37 

0.40 

0.42 

0.44 

0.46 

(last cell) 

0.31 

0.31 

0.36 

0.41 

0.40 

0.42 

0.44 

0.46 

0.47 

Axial transit-time factor 










(first cell) 

0.64 

0.86 

0.82 

0.75 

0.73 

0.68 

0.64 

0.61 

0.58 

(last cell) 

0.81 

0.81 

0.75 

0.69 

0.60 

0.65 

0.61 

0.58 

0.55 

Effective shunt impedance 










(first cell) (Mfl/m) 

27.0 

53.5 

44.6 

35.0 

39.6 

24.8 

21.5 

18.9 

16.7 

(last cell) (Mfl/m) 

47.97 

44.8 

35.2 

28.5 

25.0 

21.7 

19.0 

16. 8 

14.9 

Drift space following cavity (m) 

0.22 

0.6 

0.75 

1.0 

1.0 

1.0 

1.0 

1.0 

- 

Number of full drift tube 

55 

59 

34 

28 

23 

21 

20 

19 

18 

Average axial field (MV/m) 

1.60 

2.0 

2.6 

2.6 

2.58 

2.56 

2.56 

2.56 

2.56 


2.31 









Average gap field 










(first cell) (MV/m) 

7.62 

10.0 

8.7 

7.03 

8.9 

6.4 

6.1 

5.8 

5.6 

(last cell) (MV/m) 

7.45 

6.45 

7.2 

6.3 

6.4 

6.1 

5.8 

5.6 

5.4 

Peak surface field 










(first cell) (MV/m) 

8.9 

12.6 

13.1 

12.9 

14.0 

14.1 

14.2 

14.3 

14.5 

(last cell) (MV/m) 

10.2 

9.7 

12.9 

13.2 

14.1 

14.2 

14.3 

14.5 

14.8 

Cavity excitation power (MW) 

0.61 

1.38 

2.245 

2.48 

2.49 

2.33 

2.65* 

2.7* 

2.75* 

Total power per cavity for 100 mA (MW) 

1.58 

4.09 

5.11 

5.12 

4.88 

4.58 

4.81 

4.75 

4.68 


Total accumulated length 
Total number of unit eelb 
Total number of full drift tubes 
Total cavity excitation power 
— Total linac rf power for 100 mA 

“Estimated 


144.8 m 
286.0 
277.0 

19.6 MW 

39.6 MW 


200.30 MeV 
2.7 X 10 -3 

100 mA 

1.5 — 3.0 mrad-cro 

100 psee 
15 pps 
201.25 MHz 
400 psec 
.006 
-32° 


12 



(3.2) 



Fig. 4. Pill-box resonator accelerating cells. 

If the proton has an average velocity v and passes through the center of the cavity 
at t = 0, its coordinate is given by 


Zp — vt (3.4) 

and its kinetic energy gain while passing through the cavity is 

[9l 2 eV 0 u)Zp 
A Ek = / —-cos— -dz p 

K J-g/ 2 g v p 


sinug/2v 
= eV ° ug/2v 

(3.5) 

= eVoTtr 

(3.6) 


where the familiar sinar/x term is called the transit time factor 7f r . 
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If the proton passes through a series of cavities, each phased so that the field is 
maximum when the particle goes through the center, then cumulative acceleration can 
take place. We see, however, that if each resonator gap g is equal to /?Xo/2 where Xo is 
the free space wavelength and /3 — v/c, then 7f r is only equal to 0.64 [(sin 7r/2)/(7r/2)], 
appreciably less than unity. To improve upon this situation, for a given Vo or power per 
unit length, it is clearly advantageous to reduce the gap length g : this leads naturally 
to the configurations shown in Fig. 5. 


h-g- 


-□IK. 


□ K 


CAVITY NOSE CONES 


IK 


‘ALVAREZ DRIFT TUBE 


Fig. 5. Structures which increase the transit time factor 7} r by decreasing g. 

The cavity nose cone and the length of the gaps between Alvarez drift tubes can 
both be adjusted to bring 7} r closer to unity. Ultimately, reducing g leads to sparking 
due to excessive field gradients. There is little to be gained in T< r by reducing g to 
less than, say f3\ o/4. Table 1 shows how the gaps between successive drift tubes in 
the nine FNAL Alvarez tanks have been adjusted to handle these problems. 

3.3 Shunt Impedance 

The next question that arises is: How much power do we need to obtain a given 
amount of energy gain? Power is consumed in basically four areas: the cavity walls, 
the beam, the transmission line between the source and the cavity, and in reflections 
which decrease the power delivered to the cavity. For now, let us neglect all but the 
first. Unless one uses a superconducting cavity, appreciable power is dissipated in the 
walls of a cavity to maintain high fields in it. The power lost in one cell, Pi os t, is 
proportional to Vq. The shunt impedance R of the cavity is defined as: 

R = p 

*lost 

If the cavity is of length d, then for a multi-cavity accelerator it is also customary to 
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define a shunt impedance per unit length r where 


i PloJd 

For a negligible beam intensity and a lossless matched line between the power source 
and the cavity, one has by definition: 



- - Psource — Plost 

and the achievable value of Vq becomes 

Fq — \JR Psource 


(3.9) 


(3.10) 


where the peak energy gain, from (3.6) is 

±E K = ejRTf r Psource (3.11) 

Referring again to Table 1 for a numerical example, it can be seen that at the injection 
end of the FNAL linac, the first tank containing 55 drift tubes has an average shunt 
impedance 

R = 38 X 10 6 X 7.4 = 281.2 Megohms 


and 


T tr = 0.75 . 

Then with P SO urce — 0.6 MW 

A E k = e ^281.2 X (0.75) 2 X 0.6 = 9.74 MeV . 


3.4 Cavity Q, r/Q and Filling Time 

To calculate Pi ost it is necessary to know the quality factor Q of the cavity which 
is defined by 


q = uw *t 
Plost 


(3.12) 
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where w s t is the maximum stored energy and it is assumed that the only net flow of 
power into the cavity exactly replenishes Pi 08 t • Combining (3.7) and (3.12) one obtains 
the important ratio 


R_ v$ 
Q OJWgt 

or per unit length, 


(3.13) 


r _ (Vp/rf) 2 

Q w(ui (l /d) 


(3.14) 


Notice that, at a given frequency, this ratio is purely a property of the cavity geometry 
and does not depend on the wall material or condition, or on the quality of welds, joints, 
etc. It can be measured using dielectric beads on threads or rods. On the other hand, 
the separate quantities R, r and Q depend on all these factors. 


We will return later to a discussion of power flow along a chain of resonators, but 
for now let us notice that the natural damping of a free excited single cavity, by (3.12), 
is 


dWgt U! 

ir = -Q w “ 


(3.15) 


As a result, the time for the electric field to decay to 1/e of its initial value is 



(3.16) 


3.5 Phase Stability and Adiabatic Damping 

When a bunch of protons is injected into a linac, it inevitably has a spread of 
velocities. Superficially this suggests that in a sufficiently long accelerator, all protons 
will eventually drop out of step with the accelerating fields. Fortunately, phase stability 
prevents this from occurring. 

The important idea is that a particle whose velocity is somewhat too low gradually 
drops back in the bunch, thereby passing the structure at a phase where the acceler¬ 
ation is greater, permitting it eventually to catch up and even overtake the center of 
the bunch. It is presupposed that the drift tubes, gaps and fields are so arranged that 
one can speak of a reference or synchronous particle which is accelerated through the 
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structure without any longitudinal oscillation and defines the bunch center. In prac¬ 
tice, this is achieved by empirically adjusting the phases of the drive to the successive 
individual sections of the linac. 

It turns out to be simpler to use z, the distance along the accelerator as the 
independent variable rather than time. In this section we will use capital letters to 
represent the dependent variables for a specific particle. This will permit us to use 
lower-case letters to represent deviations away from the variables of the reference 
particle. For example, let the particle pass point z at time T. At this time the phase 
of the RF accelerating field is $ = a/TV The synchronous particle passes the same 
point at time T s and phase 4> 5 = uT s . The time and phase differences are 


/ —T - T g 

<f> =$ — = oj{T — T s ) = ut . 


(3.17) 


Referring back to Eq. (3.3), the rate of increase of energy of the synchronous 
particle is given by 

d J^- = q E !0 cos*„ 
dz 

where E z q represents the average electric field in the vicinity of z, including the transit 
time correction. For maximum acceleration, should be adjusted to zero but this 
would not yield phase stability as we shall see. For the general particle the rate of 
energy gain is 


dE K 

dz 


— q E z0 cos$ . 


The energy offset e = Er — Ej^ s satisfies 


de 

e = = q E z o [cos($ a -I- <t>) — cos <I> S ] 

~ -q Ex <j> sin <I> 5 


(3.18) 


(3.19a) 

(3.196) 


where the dot notation is intended to remind the reader that z enters the equations 
the way time does in regular mechanics. Note that we are using q for charge to free 
up the symbol e for energy offset. 
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Our intention here (and elsewhere in these lectures) is not to give the most general 
formulation but rather to explain the ideas using formulae that are as simple as possi¬ 
ble. Hence we will emphasize the linearized equation (3.19b) rather than the general, 
nonlinear equation (3.19a). For the same reason we will use non-relativistic mechanics 
even though the full relativistic treatment is straightforward (see the article by H. G. 
Hereward in LA). As it happens, for the critical early stages of a proton accelerator 
where these ideas are most important, the non-relativistic approximation is very good. 


Equation (3.19b) shows that there is an energy correction proportional to the phase 
discrepancy <j> but this is only half the story. There is also a phase correction due to 
the energy offset e since a more energetic particle has a greater velocity which allows 
it to gain in phase. T and T e are given in terms of the velocities V and V 8 by 


T = 


[ z dz ' . t - f Z dzt 

Jo V ’ 8 Jo V 8 ' 


(3.20) 


From these and (3.17) we get 



(3.21) 


It was the simplicity of this relation and the natural occurrence of 2 in (3.18) which 
motivated the choice of 2 as an independent variable. 

Velocities and energies are related by 


Ek = \™V* -, E K ,=\mVi 


(3.22) 


and as a result 


e = ±m(V*-V?) (3.23a) 

~mV,[V-V t ) . (3.23 b) 


In practice, the energy spreads tend to be very small and the linearized form (3.23b) is 
a satisfactory approximation (more so, for example, than the linearized form (3.19b) 
which won’t quite do in some cases). Combining (3.21) and (3.23b) we obtain 

= • (3-24) 
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We obtain “Newton’s law” by differentiating (3.24) and substituting from (3.21) 

oj 

=-T739^50[ c °s($a + <t>) ~ cos$ 5 ] (3.25a) 

mV a 

^ — ^qE z0 (f> sin <I> a . _ (3.256) 

m 'S 

Even the non-linear form (3.25a) has the simple feature that the “force”, “F”, depends 
only on the displacement <f>: this motivates one to introduce an “effective potential 
energy function” 


“ y ” 


q E z o [sin($ a + <{>) - <t> cos $*] . 

mb 


(3.26) 


The effective force is plotted in Fig. 6 a and the effective potential is plotted in Fig. 6 b. 
Clearly Eq. (3.25b) represents an oscillatory system provided that sin<I> s is negative. 
The “frequency” of small oscillations is given by 


2 tt_ uqE 0 \ sin $a | 

\ z \ mV? 


(3.27) 


Here, since the independent variable is z, it is the wave number 27T /\ z which ap¬ 
pears rather than u. \ z is the distance along the accelerator for one cycle of longitu¬ 
dinal oscillation to occur. Note that \ z becomes large as V s increases. Relativistically 
this effect is even more pronounced. The correct relativistic expression differs from 
(3.27) by the replacement V s —► 7 S V 3 where is the usual relativistic factor ( 1 —/? 2 )~ 2 . 

As usual with oscillators, it is enlightening to write the equations of motion in 
Hamiltonian form and to discuss the motion in phase space for which the axes are <f> 
and e as in Fig. 6 c. The Hamiltonian is given by 

__ H = -q Eq [sin ($ 3 + (f>) - <f> cos $ a ] - - 2 ~ys e 2 (3.28 a) 

9 
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Fig. 6. Phase stability graphs. 


and the equations of motion (3.19) and (3.24) are 


dll 

d<f> 


(3.29) 



(3.30) 
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For small oscillations, and dropping a constant term, H reduces to 


H = -p EzO Isin^l^-^^c 2 . (3.286) 

When a Hamiltonian does not depend on the independent variable, it is a constant of 
the motion. In that case a phase space trajectory is an ellipse with the equation 

H((j>, e) = Hi (a constant) . (3.31) 

Just as with a pendulum for which the maximum velocity can be calculated if the 
maximum height is known, or vice versa, we can obtain 

Hl~ ~ 2*1 I s ^ n $s\<i>max = ~ 2 m y3 e "iai • (3.32) 

3 

In other words, the aspect ratio of the phase space ellipse is determined by the Hamil¬ 
tonian 


e max 

fimax 



(3.33) 


But H is not independent of z. The very fact that we are discussing an accelerator 
implies that V 3 varies with z and the ideas presented here cannot be applied without 
further discussion. If the rate of variation is slow enough, the motion in phase space 
will still follow an almost elliptical path close to one of those shown in Fig. 6c. Slow 
(or adiabatic) in this context means that the fractional variation in V 3 should be small 
as the particle travels a distance \ z . It is useful to evaluate \ z for the data describing 
the FNAL linac and given in Table 1. You should find that the adiabatic condition is 
quite well satisfied. 

But in such an adiabatically varying situation, the value of e) does not remain 
even approximately constant. On the other hand, it is not hard to show (e.g., Landau 
and Lifshitz, Mechanics) that the area of the phase space ellipse is an “adiabatic 
invariant.” In other words, 


Cmax^max — constant 


(3.34) 
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Taken together equations (3.33) and (3.34), by fixing both the product and ratio of 
e-max and fimax, fix them individually. It follows that their variation is given by 


emax& 


\EzQ 


sin V; 


1/4 


(3.35) 


<t> 


max 


OC 


E z0 |sin$ s 



(3.36) 


The latter relation shows that the bunch naturally shrinks down and becomes short 
compared to the phase acceptance of the linac sections. Even the energy spectrum 
is, in a manner of speaking, damped since one is normally interested in the fractional 
energy spread e ma x/E k3 . 


We have not said much about the motion for large amplitudes where (3.28a) must 
be used. Nor will we, as it is assumed that the reader has already encountered phase 
space plots like Fig. 6c elsewhere. On injection, particles inside the separatrix are 
captured. Those outside are largely lost although even some of them can be captured 
because of adiabatic damping. 

3.6 Transverse Defocusing 


There is one undesirable consequence of the phase stability considerations of the 
last section. We saw, following Eq. (3.25b), that longitudinal stability requires that 
sin be negative. Referring back to Eqs. (3.2) and (3.3), the time rate of change of 
the electric field is 

dE Z jr, 

—— = — E z o oj cos u>t 

at 

which, for the synchronous particle becomes 

^ = -E z0 ojs'm^ 3 . (3.37) 


This expression is positive if there is longitudinal stability. Looking at Fig. 7 it can 
be seen that fringe fields in the gap region lead to transverse forces. As the proton 
enters the gap it feels a focusing force and, as it leaves, a defocusing force. There is 
a net defocusing effect since from (3.37) the magnitude of the field is increasing with 
time. 
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A quick way of seeing the inevitability of this effect is to consider the motion in a 
moving frame of reference in which the synchronous particle is at rest. In this frame 
there is only a time-independent electric field since the phase velocity of the field is 
equal to the particle velocity. In this frame, Earnshaw’s theorem which applies to DC 
fields and potentials excludes the possibility of a potential minimum, and motion in at 
least one of the three dimensions must be unstable. Since the longitudinal motion has 
been adjusted to be stable, the transverse motion will not be. 



Fig. 7. Transverse focusing fields in accelerating gaps. 


It is possible to show this effect mathematically. So far we have not explicitly 
given the 2 -dependence of the accelerating field other than to say that the protons are 
inside drift tubes where there is no field except at times when the field is in the right 
direction. As mentioned above, we can represent this by a traveling wave given by 

E z (z,t) =E z0 cos . (3.38) 

The argument has been arranged so that a particle of speed V s stays at a fixed phase 
of the traveling wave. In Sec. 4 we will see how, by Fourier transformation, such a 
traveling wave can be obtained for a general accelerating structure, but for now let us 
accept (3.38) as the dependence on z and t separately, for example near the center of 
one of the accelerating gaps. 

We will now use this expression of the field to find the form of the radial Lorentz 
force equation. Taking r as a transverse radial displacement and 6 as an azimuthal 
angle, we can show that, for small r, the form of (3.38) and Maxwell’s equations fix the 
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other field components. We also assume that H z vanishes and, assuming rotational 
symmetry, that Eg = H r — 0. In the absence of free charge, Gauss’s law makes V -E 
vanish, i.e., 


r dr * E ) dz 


(3.39) 


For small r, E r will be proportional to r, and we get 

^=5 £ (rft) = - £ E ’°V, sinu { t ~Io^) ■ 


(3.40) 


By Ampere’s law 


and from this we obtain 


1 d_ 

r dr 


(rH e ) = Cq 


dE z 

dt 


TW *4 

Hf = /■.,() ftisin 0 J 



(3.41) 


(3.42) 


We can now write the radial Lorentz force equation for the synchronous particle 


A( m f) = q(E r - VsHoHg) 


quE z o 

2V 8 


(l-^f)rsin $ 3 


(3.43) 


where it has been recognized that the phase factors in (3.38), (3.40) and (3.42) are 
the same as in (3.25b), i.e., But since the sign in (3.43) is opposite to that in 
(3.25b), one or the other must lead to unstable motion which is what we set out to 
demonstrate. It should be noticed though, by the presence of the factor 1 — , that 

the defocusing vanishes in the relativistic limit. 

Space charge forces also tend to enlarge the transverse beam size but the effect 
just described limits operation even at low beam current. Early corrective schemes 
using electrostatic grids were not very satisfactory. It was the discovery of alternating- 
gradient focusing by Blewett, Courant, Livingston and Snyder and the introduction 
of quadrupole magnets inside the drift tubes which solved the problem. 
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3.7 Applications and Some Practical Aspects 

For high energy physics, the main application of proton and ion linacs is as injectors 
into circular machines, but by themselves they also have numerous applications in 
nuclear research and medicine. Included are meson factories (LAMPF), heavy ion 
research (SUPER HILAC, UNILAC, RILAC) and intense neutron sources. 

A typical layout of a proton linac is shown in Fig v 8._ Th^source is generally a 
metal envelope containing a gas (hydrogen or other) which is ionized into a plasma. 
The protons (or ions) are formed into a jet through a cylindrical nozzle. 


200MHz 200 MHz 


800 MHz 


DPHC W 


—[bTHb^I— | ALVAREZ | ---|SJDE : Cp_UP_LED_S.j 

■ ■ 1 ■ 


750 kV 


/ V 

MAI 


FNAL 

BNL 


LAMPF' 


Fig. 8. Typical proton or ion linac layout. 

This system, called a duoplasmatron (DP), floats at a large negative voltage supplied 
by a Cockcroft-Walton supply and emits the protons (or ions) into a high-gradient 
accelerator column (CW). This results in particles with an energy of as much as 1 
MeV (typically 750 keV) but near ground potential. Photographs of the BNL versions 
of these components are shown in Figs. 9 and 10. 

This pulse of protons forms a DC current of perhaps 100 or 200 mA for 100-200 
microseconds. The protons then become bunched in the drift space following the 
low power, velocity-modulating, 200 MHz buncher labelled B\. The second buncher, 
£?2> used for longitudinal phase space matching, injects the protons into the Alvarez 
linac where the capture efficiency at a synchronous phase = —30° is as much as 
80%. Figures 11, 12 and 13 show successively more magnified views of the BNL linac. 
After 150 or 200 m the protons have been accelerated to about 200 MeV with a pulse 
repetition rate of 10 or 15 pps. The parameters of other injectors, such as that at 
FNAL, are similar. At LAMPF this stage is less than half as long and achieves 100 
MeV with a greater repetition rate (120 pps) and a 500 //sec pulse length. 
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4664A9 

Fig. 9. Cockcroft-Walton pre-injector (750 kV) 
(BNL-HAEFELE) 
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Fig. 10. High-gradient accelerating column (750 kV). 

(BNL) 







Fig. 11. BNL Alvarez linac from low-energy end 
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Fig. 12. BNL Alvarez linac showing drift tubes, multi-stem 
arrangement, tuning bars, and slugs for field stabilization. 
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Fig. 13. First Alvarez drift tubes (notice sparking). 

In all these Alvarez linacs, typical average axial fields are in the range of 1 to 2.5 
MV/m, with gap fields at 6 to 10 MV/m. There are many tanks, each containing 
many drift tubes (refer again to, e.g., Table 1). Each tank is supplied with a peak 
power of several MW, typically ~ 5 MW. 

To obtain higher efficiency, the remainder of the LAMPF accelerator which goes 
up to 800 MeV utilizes a structure at higher frequency, 800 MHz. This structure uses 
side-coupled cavities. We will defer explanation of this important development until a 
more systematic formalism for periodic structures has been developed in Sec. 4. 

3.8 Recent Developments 

One recent invention, first proposed in 1970 by I. M. Kapchinskii and V. A. 
Teplyakov at ITEP in the Soviet Union, is to replace the front end of the linac struc¬ 
ture by a device in which the RF fields that are used for bunching and acceleration 
can also serve for transverse focusing. In the USA, the structure which can perform 
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these combined functions is called the RFQ (Q for quadrupole). As mentioned in Sec. 
2, RFQ’s have now operated successfully at several laboratories. Table 2 gives param¬ 
eters for the designs of three such devices. 


Table 2. Design Parameters for Three Different RFQ’s 

Los Alamos CERN BNL2 

Polarized 


Ion 

H~ 

H+ 

//- 

Frequency (MHz) 

425 

202.56 

201.25 

Input energy (MeV) 

0.10 

0.05 

0.02 

Output energy (MeV) 

2.00 

0.52 

0.75 

Nominal current (mA) 

100 

100 

0.1 

Transmission efficiency (% ) 

91.4 

93.6 

100 

Output emittance a 

0.021 

0.05 

0.011 

Length (cm) 

288.6 

138.2 

123.4 

Peak rf power (kW) 6 

1020 

304 

66 


°RMS normalized area/jr in cm mrad units. 

^Includes power required for the rf manifold and for the beam. 


To understand how the RFQ works, refer to the early model built at ITEP (Fig. 14) 
and the simplified sketch of the vanes which make up the four opposing arrays of poles 
(Fig. 15). 

The fields and currents are said to resonate in a modified TE 210 mode. At any 
particular distance along the resonant structure, the poles are shaped and powered to 
give an electric quadrupole field. As the fields vary sinusoidally in time, an alternating 
gradient with net transverse focusing results. However, since the poles are also given 
a sinusoidal-like radial variation, this configuration leads to longitudinal acceleration 
and adiabatic bunching as well. 

One reason for which this device is of such great interest is that in a small space, 
it can replace the very large Cockroft-Walton section of an ordinary injector, as can 
be inferred by noting the typical input and output energies in Table 2. Although the 
psak RF powers are relatively high, the RFQ’s hold the promise of far better capture 
efficiency and reduced emittance growth at reduced cost. 
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Fig. 15. Transverse and longitudinal views of the poles of an RFQ 


Fig. 14. RFQ structure built at ITEP, Moscow. 


4664A14 
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c-c 


Schematic drawing of RFQ resonator. 



FMIT Prototype RFQ 


Ion: Deuteron 



- j WM 

□HUE 

'—- GB 

-HUE 


Frequency: 80 MHz 

W(MeV) 

0.075 

0.075 

0.120 

0.80 

2.00 

Nominal Current 

: 100 mA 

E (MV/m) 

8 

3.01 

17.6 

17.6 

17.6 

17.6 

Current Limit: 

205 mA 

a 

1.00 

1.00 

1.12 

2.13 

2.13 

PARMTEQ RESULTS 


a (cm) 

8.31 

1.42 

1.35 

0.89 

1.00 

Input Current: 

106 mA 

r o (cm) 

8.31 

1.42 

1.42 

1.42 

1.63 

Output Current: 

100 mA 

$ s (deg) 

- 

-90 

-72 

-30 

-30 

raj Input (90!t): 

0.070 

E (MV/m) 
o 

0 

0 

0.40 

1.97 

1.31 

raj Output (905i): 
*-* 

0.160 

V (kV) 

185 

185 

185 

185 

185 

■g| Input (RMS): 

0.017 

B 

0.20 

6.82 

6.82 

6.82 

5.20 

Output (RMS): 

0.036 

L(m) 

0 

6.7 

130 

274 

388 


CELL NUMBER 



FMIT prototype RFQ pole-tip shape. 

Fig. 16. Characteristics of RFQ for FMIT linac (LANL). 
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Figure 16 gives detailed design information for the RFQ for the LANL FMIT 
project to accelerate deuterons. The cavity containing the 3.88 m long vane shows the 
pole profile in the four successive sections of the device: the radial matching section 
(RM), the bunch shaper (S), the gentle buncher (GB) and the accelerator (A). 

Another interesting modern development is the linear induction accelerator men¬ 
tioned in Sec. 2. The principle of this accelerator is illustrated by the module shown 
in Fig. 17. 


LINEAR INDUCTION ACCELERATOR 



As stated earlier, the idea of the induction linac is based on Faraday’s law. The 
accelerating action can be described as analogous to a transformer. The primary of the 
transformer is excited by a high voltage source connected to a pulse forming network, 
switch and transmission line. The increasing current flowing through the primary 
when the switch is fired creates a rapidly changing magnetic flux in the toroidal core 
(thin-laminated iron or ferrite) which in turn induces a voltage across the accelerating 
gap. The beam then can be seen as the secondary winding. The accelerator itself 
consists of an array of such modules which can be triggered in succession at a rate 
commensurate with the velocity of the particles, their particular q/m and the mean 
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accelerating field. Accelerating capability, expressed as the product of gap voltage 
times pulse length, is limited by the saturation of the ferromagnetic material, which 
in turn is determined by the amount and type of material present. Typical average 
gradients of 2.5 MV/m can be reached. The technology of the induction linac seems 
to be well suited for very high beam currents (100 to 10,000 A) with short pulses (10 
to 2000 nsec) at repetition rates up to 100 Hz. 

- The first version of the induction linac was the Astron injector designed and com¬ 
pleted by N. Christofilos at LLNL in 1963. More advanced versions were built in 1971 
at LBL and Dubna for their respective electron-ring accelerator (ERA) projects. The 
idea of the ERA was to form and accelerate rings of electrons and to embed protons 
or ions in the potential well of the electrons circulating in the rings. By accelerating 
the rings to relativistic velocities without “spilling” the protons, it was hoped that 
extremely high-gradient proton machines could be built. Unfortunately, it was discov¬ 
ered after some time that instabilities would prevent this scheme from being realizable 
at high energy. As a result, the ERA projects lost their momentum although at Dubna 
this type of accelerator is still used at low energy for ions. In the process, however, the 
feasibility of practical induction linacs for very high currents was established. Table 3 
gives the parameters of some induction linacs in existence. All of them happen to be 
for electrons. 

In recent years, the induction linac has been proposed for a new application: heavy 
ion fusion. In this application, the linac which is called the driver is to produce pulses 
of ions such as, for example, Krypton or Thallium (see Fig. 1) on the order of three 
megajoules. This pulsed beam energy, to be achieved by accelerating the ions up to a 
maximum of 20 GeV and currents of 15,000 A, is to be focused simultaneously from 
several directions onto a deuterium-tritium pellet to produce what is called inertial 
confinement fusion. The pulse length of the induction linac modules can be reduced 
along the machine as the energy of the ions increases and the time of passage through 
them decreases. Typically, the pulse length might be 2000 nsec at the front end and 
20 nsec at the output. Progressive beam bunching is obtained by shaping the current 
pulses from the pulse forming networks so as to produce a positive voltage ramp 
within each module. Thus, later particles with higher energies catch up with earlier 
ones at lower energies. Even at an output energy of 20 GeV, the particles are still 
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non-relativistic and final power amplification up to several tera-watts can be obtained 
by longitudinal bunch compression in the transport lines leading to impact on the 
pellet. 


Table 3. Parameters for Typical Induction Accelerators 


Accelerator 

Astron Injector 
Livermore 
1963 

ERA Injector 
Berkeley 
1971 

NEP 2 Injector 
Dubna 
- 1971 " 

ATA 

Livermore 

1983 

Kinetic energy, 
MeV 

3.7 

4.0 

30 

50 

Beam current on 
target, A 

350 

900 

250 

10,000 

Pulse duration, 

300 

2-45 

500 

50 

ns 





Pulse energy, 
kJ 

0.4 

0.1 

3.8 

25 

Rep rate, pps 

0-60 

0-5 

50 

5 

Number of 

300 

17 

750 

200 


switch modules 
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4. Particle Acceleration by Guided Waves 


We will now present a somewhat more general discussion of particle acceleration 
by guided waves. As we have seen in the previous section, a series of appropriately 
phased accelerating gaps, separated by drift spaces, can be analyzed simply to give 
a satisfactory description of a proton linac. The present discussion will include such 
structures as a special case but will focus on the guided wave -itself rather than on 
the details of the structure. In the interest of simplicity we will concentrate mainly 
on particles with speeds near the velocity of light (i.e., electrons) although the same 
formalism applies also to slower particles. Rather than building on our previous dis¬ 
cussion of proton accelerators, we will begin again and consider the general problem 
of particle acceleration by an electromagnetic wave. 

4.1 Motivation 

As mentioned before, we are concentrating on particles with velocity near c and this 
makes it natural to attempt to use an electromagnetic wave in free space since it travels 
at the same speed and might have a cumulative effect on the particles. Unfortunately, 
such a wave has the property that its E and //-field are orthogonal in space and have 
components only in the transverse plane perpendicular to the direction of propagation. 
Since, in addition, the wave is bound to be of limited transverse extent, only one of 
two outcomes is possible. Either the particles start along a direction exactly collinear 
with the wave, in which case they see no force at all (i.e., the transverse electric and 
magnetic forces cancel exactly); or they start at a slight angle, in which case the 
transverse component of the E'-field produces some acceleration along the direction of 
motion of the particles but the effect is not cumulative because the particles soon find 
themselves outside of the field of the wave. 

There are two ways around this difficulty: either one keeps bending the beam back 
into the wave or one reflects the wave back into the beam. The first is the principle of 
the Inverse Free-Electron Laser which is discussed elsewhere in these Proceedings. The 
second one is the principle of all linacs using rf waveguides. An appropriate waveguide 
can redirect the electromagnetic wave to keep it superimposed on the particle trajec¬ 
tory. At the same time, it can also be such that it has a longitudinal electric field in 
the direction needed to impart energy to the particles. To have a cumulative effect, 
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however, the particles and the wave must be synchronous. As we will now see, it is 
the need for the simultaneity of these conditions which presents the main challenge to 
the accelerator structure designer. 

We will now go through a short review of guided waves to clarify these points. 

4.2 Wave Confined by Parallel Planes 

Figure 18(a) shows two plane waves propagating inTree space at angle 9 on either 
side of the 2-axis, and with their F-fields lying in the xz plane (the plane of the paper). 
The lines represent a snapshot of wavefronts of maximum | E | at a fixed time. They 
are labelled +, —, +, etc., to indicate whether ft points left or right. The free space 
wavelength Xq is shown. It is given by 

X 0 = c— . (4.1) 

w 


At a few of the intersections of the wavefronts the vector sums of the fields due to 
the two waves are shown. On the dashed lines the JS’-field is transverse to the 2 -axis. 
Along these lines it is permissible to place a mirror (i.e., a conducting plane) without 
affecting the fields. 

Now the two waves can be thought of as bouncing back and forth between the two 
conducting planes. Changing our perception somewhat, we can see in Fig. 18(b) that 
there is a single wave having a longitudinal component of & and propagating in the 
direction of the 2 -axis. The guide wavelength \ g is shown. It is the distance along the 
2 -axis after which the pattern repeats and it is given by 

, _ ^0 
9 cos 9 



as can be deduced from the geometry of the figure. 

We are now confronted by one of the obstacles which tends to foil our efforts. If we 
wait for a time 2n/ui which is one period of oscillation, the pattern will have advanced 
a distance \g. That is, the phase velocity is 


u c 

= X„ — = X„ — 


2?r 


^0 


c 

cos 9 


(4.3) 
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DIRECTION OF THE 
TWO PLANE WAVES 



Fig-_18. (a) The origin of longitudinal electric fields in guided waves, and (b) their 

configuration between two parallel conducting planes. 



Since cos 9 is less than 1, the phase velocity is greater than c. Since the particle 
to be accelerated is necessarily traveling at a velocity lower than c, our wave cannot 
yet be used for cumulative acceleration. 

Before addressing this problem, we should discuss some other preliminary questions 
regarding waveguides. The conducting wall could be moved out to any of the other 
planes on which E z vanishes and this introduces the possibility of different waveguide 
modes . Also, it is natural to “close” our waveguide, for example by making it rectan¬ 
gular through the addition of conducting walls parallel to the plane of the paper: this 
makes even more modes possible. In this case it is still possible, although harder, to 
visualize our wave as made up of plane waves reflecting off the walls - four waves must 
be used. Rather than attempting this, we shall withdraw to a more abstract treatment 
based on Maxwell’s equations. The wave we are considering has a transverse E x and 
a longitudinal E z component but only a transverse H y magnetic component, which 
makes it a so-called TM wave. All field components can be derived from E z which 
satisfies the wave equation 


d 2 Cz , 
dx 2 


or 


^E z = -- I 


E, 


(4.4) 


where we have already built in some assumptions which we will now spell out. We 
have assumed a time dependence of the form eP* 1 and a ^-dependence e~^ z where 7 is 
complex: 


7 — a + j/3 . (4.5) 

Also, for simplicity, we are limiting ourselves to the case of Fig. 18 with no y- 
dependence. By inspection, we guess a solution satisfying (4.4) and the boundary 
conditions at x = 0 and x — a, with amplitude C, 


E z 



(4.6) 


For now, a is the guide width. Substituting (4.6) into (4.4) we obtain, for this simplest 
possible case, the relation between guide wavelength and frequency 



(4.7) 
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For other modes the first term is changed, but for simplicity we will not generalize 
the notation to cover all cases. The reader should be aware of the fact that some of 
the following formulae, though correct in form, may need appropriate modification for 
other modes. 

If there is no attenuation, a vanishes and (4.7) yields 



where ft is the wave number corresponding to the guide wavelength \ g \ 



(4.8) 


(4.9) 


This ft (which should really be called ft g or k g ) must not be confused with the usual 
ft — ti/c used for relativistic particles. Expression (4.8) is called a dispersion relation. 
It can be rewritten in terms of the “cut-off frequency,” a term to be justified shortly, 
given by 

oj c = — c (4-10) 

a 

whereby it becomes 

c 

The graph of u versus ft in this field is called a Brillouin diagram. It is easily recognized 
to be a hyperbola (see Fig. 19). When w drops below u c , (3 becomes imaginary and 
there is no more propagation. That is why u c is called the cut-off frequency. 

Following the usual analysis of a dispersion relation, the phase velocity is given by 


_ u> c 

P P \A “ (We/w ) 2 


(4.12) 


which, as expected, is greater than c, as can also be noted from the slopes shown in 
Fig. 19. 
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Fig. 19. Brillouin diagram for guided waves in an unloaded guide. 

It is useful to think about the situation at the extremes of this graph. At high u, 
(4.11) becomes the free space linear relation between frequency and wave number. In 
Fig. 18 this corresponds to wavefronts that are very close together (\ g —>■ \q, 0 —* u>/c ) 
and the wave travels freely down the guide, like a flashlight beam in a hollow tube. 
We have already seen that, although the speed v p is close to desired, the wave is not 
suitable for acceleration because its i^-field is transverse. 

At the low frequency extreme, the lines in Fig. 18(a) are as wide apart as they can 
be while still allowing the picture to be completed; that is, the spacing from + to — 
should be a, 



which agrees with (4.10). The wave number 0 is very small, corresponding to large 
guide wavelength and large phase velocity. In this case the wave bounces back and 
forth at right angles to the guide; it is not useful for acceleration because of the high 
value of v p . 

Since accelerator guides are generally circular, perhaps more time than is justified 
has already been devoted to rectangular guides and parallel planes in particular. But 
there are a few more points which can illustrate our ideas in this simple case, free from 
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the intrusion of complicated geometry and mathematics. These have to do with power 
flow and dissipation in the walls. 

The non-vanishing transverse field components between infinite parallel planes cor¬ 
responding to E z in (4.6) are given, using Maxwell’s equations, by 


(4.13) 

* dx 7r v a ' 

The power flow P( r is obtained by integrating the average Poynting vector over a 
transverse plane 

P tr =\ Re f E x Hydxdy . (4.14) 

Since we are assuming that there is no y-dependence, i.e., the parallel planes are infinite 
in y, we don’t have to satisfy the boundary conditions on any bottom or top plane. 
For a height b, with 7 = jft taken as pure imaginary, we get 

Ptr = C 2 ^ eou>a 2 fl . (4.15) 

The standard method for calculating power loss in the walls (an issue of dominant 
importance in linear accelerators) is to take the fields obtained above for the lossless 
case to calculate the wall currents. In our case, the loss is in the side walls where, by 
Ampere’s law, the longitudinal current I z per height b is given by 


E x = - 


7 a 


dE z 


7 a 


f**\ - 


7T 


2 dx 


= - —<7 cos — e 
7 T v a 1 
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I Z = Hyb . 


Although the model is not conceptually quite correct, one obtains the correct answer 
if one assumes that this current flows uniformly in a skin depth 6 given by 

( 2 \V 2 

« = (-) ( 4 - 16 ) 

where fio = 4x X 10~ 7 henry/m and <7 is the conductivity in mhos/m. The average 
power lost in a short length dz of height 6 of the two side walls is of the form 2 X 2 |/| 2 P, 
i.e., given by 

dp,„< _ |gj,6| 2 _ 

dz 08 b atii r 2 
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where (<x<$) 1 is the skin resistance R s . To account for this loss, the propagation 
constant must have a small real part given by 


dPio 3 t/dz _ 1 4tt R s 
2 p tr P a x 0 01o/e 0 


(4.18) 


Essentially the same information can be quantified J>y defining a Q factor for the 
guide * - 


u w 8t 

V dRlcst/dz 


(4.19) 


where w s t is the stored energy per unit length. This quantity can be calculated from 
the maximum electric field, averaged over a wavelength: 


™st — 


£o 

4 



f 0 

4 



C 2 ab 


(4.20) 


which yields 

^ _ *^o[(^r) + 1] \Jv o/ f o 
Q ~ 8 RTa 


(4.21) 


Finally, there is one more relation between these quantities. If the energy travels 
with the group velocity v g , then 


Ptr - tig W st 


(4.22) 


where 

dui Ptr 
Vg ~ d/3 w s t 


(4.23) 


Note that in Fig. 19, the group velocity is the slope of the ui — (3 diagram. It is 
equal to zero at cutoff (u = u c ) and equal to c when u —*■ oo. The reader can verify 
that both definitions of v g are equivalent and that Q and a are related by the equation 
a = oj/ 2 v g Q. 
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For a practical accelerator structure, all of the above quantities would have to be 
calculated. In what follows we will not attempt this as it can become quite complicated. 
Conceptually, however, these calculations are not different from those which have just 
been done. The exact formulae, if they can be obtained in closed form, resemble the 
ones we have just derived as far as dimensional factors are concerned. The numerical 
factors, of course, depend on the exact geometry. 

4.3 Circular Waveguide 

Most linear accelerator structures have configurations with circular cross sections. 
To simplify matters, let us take the case where there is just one boundary, r = b. If 
we limit ourselves to TM modes with no ^-dependence, the wave equation of (4.4) can 
be replaced by 


1 BE, _ 2 „ a, 2 

a o H— r 7 — 

dr* r or c* 


The other non-vanishing fields are given in terms of E z by 


Ha = 


Mo dE z 


kf dr 


(4.24) 


(4.25) 


where 



dE z 

dr 


(4.26) 


*? = 1 2 + (y ) 2 ■ (4.27) 

The simplest solution of (4.24) is the so-called TMq j mode where the first subscript 
( 0 ) denotes zero-0 variation and the second ( 1 ) denotes one radial variation. For this 
mode the fields become 


E z =C Jo{k c r)e~l z 

He = 3 ~t^ C J\{k c r)e~ lz 
Kc 

E r =lc J\(k c r)e~^ z 


(4.28) 
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If one compares expressions (4.28) with those of (4.6) and (4.13), one sees that the 
main changes are that Bessel functions have replaced sine and cosine functions. The 
waveguide radius becomes involved because of the requirement that E z vanish at r = 
6 , which means that k c b must be the first zero of Jo(k c r): 


k c b = 2.405 . (4.29) 

The guide dispersion relation takes the form of (4.11) if the cutoff frequency is defined 

by 


u c = k c c — 2.405^ . (4.30) 

Note that, as claimed earlier, this expression differs from (4.10) only by a dimensionless 
numerical factor, in this case 2.405 instead of ir. All the other waveguide properties 
such as Ptr, Plan, a, Q can be calculated as for the case described in Sec. 4.3. 
We will not make the effort to do this because the TMq\ mode still suffers from the 
shortcoming that the phase velocity, as given by (4.12), is greater than c, making it 
useless for cumulative acceleration. However, having obtained the field expressions for 
the circular waveguide, we can profitably return to our cylindrical resonator originally 
discussed in Secs. 3.2, 3.3 and 3.4. 

4.4 Cylindrical Resonator 


In the circular waveguide, the dispersion relation (4.11) determines the longitudinal 
propagation constant as a function of oj: 



(4.31) 


If however the guide is capped off at both ends by conducting planes as shown in Fig. 
20 , a constraint is placed on j3, and w can only assume discrete values corresponding 
to resonances. As in Fig. 18, the boundary conditions at the plane surfaces can be 
satisfied by superimposing two waves of appropriate wavelength, one going left, one 
right. The result is a standing wave. A simple possibility is to take the mode where 
E r — 0 at the plane surfaces and 
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i.e. 


> 


Fig. 20. Cylindrical resonator. 



Using (4.11) and (4.30) this yields for the resonant frequency 



(4.32) 


This is the so-called TMq\\ mode which has one variation in the 2 -coordinate. It 
turns out, however, that in a cylindrical waveguide, there is an even lower mode (the 
one used in Alvarez linacs) in which E r vanishes not just at the ends, as required by 
the conducting planes, but everywhere. For this mode \ g is infinite, which makes /? 
zero. The resonant frequency is 


= 2.405^ (4.33) 

which is independent of d. This mode is called the TMqio mode. The expressions for 
the field given by (4.28) are reduced to 

E z =CJ 0 (k c r) 


i h(k c r) 


(4.34) 


where rj — (^oAo) 1 ^ 2 = 377 ohms. From these, using the definitions of the shunt 
impedance R (3.7) and Q (3.12) or (4.19), it is possible to show that 




(4.35) 
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and 


q = V 2.045 
V R a 2(1 + 6 /d) 


(4.36) 


where = 1 /<t6 is as defined previously. It would be instructive for the reader to 
verify these expressions and to get a feeling for the magnitudes involved by working 
out some numerical examples for R, Q, R/Q and Tp at 200 MHz (d = 8 m) and 2856 
MHz (d — 3.5cm). 


Note that for this calculation we have neglected the effect of the holes at the ends 
of the resonator, needed to let the particles pass. As we discussed earlier, a linear 
accelerator can be built by stringing many such cavities together and adjusting their 
phases to the beam to give cumulative energy gain. However, rather than assuming 
that these cavities are independently driven, we will now try to analyze the model 
where they are connected together into a periodic chain of coupled resonators. 


4.5 Wave Propagation in Coupled Resonator Chains 

For simplicity we will start with the simplest possible resonator, namely an LC 
circuit. The pillbox resonator in the TMq\q mode is indeed little more than an LC 
circuit where C a is the capacity between the end plates and L is the inductance of 
the toroidal outer region where the magnetic field is predominant. Representing the 
coupling between resonators by a parallel capacitance 2 C p , the circuit diagram for the 
chain is shown in Fig. 21 . 



Fig. 21. Chain of lightly coupled resonators. 

Note that in the limit of large C p , the resonators are uncoupled. For the pillbox 
resonators this would correspond to having no beam hole. To insure that we are close 
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to that situation, we will assume that 


» C ( . (4.37) 

We do not intend to imply that representing the coupling by a capacitance is an obvious 
model for the actual coupling, but it is simple and leads to the same qualitative results 
as more complicated models. 

It was L, Brillouin who first systematized the analysis of such periodic chains, but 
the essential results we require can be derived in a few lines. Applying Kirchoff’s 
current law to the nth loop in Fig. 21 yields 

+ jk~ + iuL ) +i - l) = 0 

or 

*n+l - 2 cos {/3d)i n + »„_! = 0 (4.38) 

where 

cos (fid) — l +~ — w 2 C p L . (4.39) 

t-'s 

From the point of view of the present discussion we have simply introduced a new 
symbol (/3d) which is a certain combination of u and the circuit parameters. But, 
naturally, our choice of symbol anticipates identifying this quantity with the previously 
defined quantities (3 and d. 

Two solutions of the difference equation (4.38) are 

sin 

i n = (n/3d + <j> o) (n = 0,1,2,...) (4.40) 

COS 

as can be checked with standard trigonometric identities. If j3 as defined by (4.39) is 
real, then (4.40) will represent standing wave solutions which the structure can support. 
Remember that it is implicit in an impedance equation such as (4.38) that the time 
dependence e ,UJt is assumed and the real part is to be taken. By superimposing the 
two solutions (4.40) with proper phase, unattenuated traveling waves going in either 
direction can be formed. The condition for f3 to be real is that 

—1 < cos /3d < 1 (4-41) 
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which can be re-expressed as a band of propagating frequencies 



WQ < u < W* . 

(4.42) 

The frequencies wq and w^ defining the limits of the band 

are obtained by combining 

(4.39) and (4.41) to give 




1 

Wo = ,-- 

y/LCg 

- - (4.43) 


1 f C s \ 

(4.44) 


- s/lcX 1 + c p ) ■ 


Approximation (4.37) was used in obtaining w*. Note that the subscripts on wo and 
u n are the corresponding values of the angle /?d at the band edges. 


The Brillouin diagram for this structure is shown in Fig. 22 which is a graph of 
(4.39). The band of propagating frequencies includes the resonant frequency 1/ \JLC S 
of the uncoupled resonators. For weak coupling ( C p much larger than (7 5 ) the band is 
narrow, and for strong coupling ( C p closer to C 3 ) it is broad. In fact, we have 


Aw Wfl- — wo C s 

wq wo C p 


(4.45) 
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Certain features of this Brillouin diagram are important properties of periodic 
structures in general. It is even in /3d and periodic in /3d with period 2 j r since these 
are both properties of the function cos/3d. The dashed extensions in Fig. 23 convey 
this effect. But this extension is superfluous for the system we are discussing since, 
in (4.40), shifting /?d by 2n does not alter the solution. This is due to the fact that 
our single cells are fully described by one variable (the loop current) and continuous 
variation of the argument is not defined. On the otherTiand, our single cell could be 
more complicated. For example, suppose it had a mode of oscillation at some higher 
frequency as well. This would lead to another propagating band as shown in Fig. 23. 
For linear circuits, these modes do not influence each other. For reasons to which we 
will return, we have chosen to emphasize the band — ir < /3d < 0 as the principal 
band for the higher mode. 



Fig. 23. Brillouin diagram for a chain of cells, each having two modes. 


As usual with dispersion relations, the Brillouin diagram can be used to extract 
phase and group velocities. Taking d as the physical length of one cell and with oj/c as 
the vertical axis, the slope of a 45° line, as shown, corresponds to a phase velocity equal 
to c, which is required for an electron accelerator. This condition fixes the operating 
frequency at the intersection point. Another general feature is that the group velocity 
(given by the slope doj/d/3) is always less than c and vanishes at ujq and at cj n . 
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4.6 Periodically Loaded Structures 

We now have to face up to the fact that a linear accelerator is more complicated 
than a chain of one-dimensional resonators. The individual cells support not just one, 
but many modes and it may not be justified to treat the coupling as weak. While the 
latter fear is somewhat unfounded, it still makes sense to attack the problem from yet 
another point of view - analyzing the structure as a periodically-loaded waveguide. 

. This-brings us back to the question which was left open in section 4.2: how do 
we slow down the phase velocity of a waveguide? For the chain of pillbox resonators 
the trick was easy since in the lowest mode the frequency is independent of the length 
d and in the weak coupling limit the resonators can be phased independently. There 
are many other ways of slowing down the phase velocity, some of which are shown in 
Fig. 24. 

Of these structures, the simplest is the iris-loaded or disk-loaded waveguide. This 
is just another name for the pillbox resonator chain which we have been discussing, 
except that now the iris hole radius, a, has become larger. Treating the irises as 
perturbations, we can ask what effect they have on the waveguide modes. If there is a 
wave propagating along the guide in the +z direction, there will be a reflection at each 
iris. For a band of frequencies, these reflections will interfere destructively because 
of the phase shifts resulting from the spatial separation of the irises. As a result, at 
these frequencies the irises will have only a minor effect on the propagation and the 
dispersion relation will resemble Fig. 19. 

But there are certain frequencies, for which the reflected waves from successive 
irises are exactly in phase. At such frequencies, unattenuated propagation is impossible 
since reflections necessarily result in comparable amplitudes for forward and backward 
going waves. Standing waves are possible, however. The condition for this “Bragg” 
reflection is 


^ = ~ n n =0,1,2,... (4.46) 

since this makes the phase advance in going ahead to the next iris and then back an 
integral multiple of 2x. This situation is again characterized by a vanishing group 
velocity. 
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We can also anticipate that for frequencies near w n7r there will be a strong effect 
of the irises on the dispersion relation between w and /?. To investigate this, consider 
Fig. 25 in which we have guessed two standing-wave field configurations of opposite 
symmetry relative to the irises. In the figures, lines and arrows indicate electric fields, 
and circles indicate magnetic fields. It can be seen that in case (b), since E r vanishes 
at the locations of the irises, the presence or absence of the irises does not have much 
effect. For this case then, the relation between u and ]3 will be” very similar to that 
for a waveguide without irises. But in case (a) the fields cannot be as drawn since the 
presence of the irises must force E r to vanish there, at least at the radial positions 
occupied by the metal irises. In other words, currents flowing in the irises cause the 
standing wave of wavelength 2d to have a different (actually lower) frequency than if 
the irises were not there. 



(a) 



Fig. 25. “Guessed” standing-wave field configurations of opposite symmetry with 
respect to the irises at the stop band. Lines and arrows indicate electric fields, circles 
indicate magnetic fields. 


We now have in hand all the elements required to complete qualitatively the Bril- 
louin diagram for the disk-loaded waveguide (Fig. 26). We can start with Fig. 19 for 
the unloaded guide. As /? increases away from zero, w will initially follow the hyperbola 
but will then flatten out and become horizontal at /? = n/d where the group velocity 
must vanish. We have just seen that there is another frequency giving this same wave 
number and that, furthermore, it lies on the hyperbola for the unloaded guide. This 
allows us to complete the branch between /3 — n/d and ft = 2n/d. 
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Fig. 26. Brillouin diagram for a disk-loaded waveguide. 
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Similar considerations lead to successively higher branches. Again, for an electron 
accelerator, possible operating points lie on the 45° line. Clearly, the loading disks are 
capable of lowering the branches to intersect this line as shown. For normal operation, 
the intersection with the lowest branch is chosen. It is indicated by uq p , fop- 

4.7 Space Harmonics 

A more quantitative understanding and completion o£Fig.-26 requires that one 
find a method to calculate the exact dispersion relation and the corresponding electric 
and magnetic field expressions for a given slow-wave structure. To match the periodic 
boundary conditions shown in Fig. 27, it is customary to invoke a theorem ascribed 
to the French mathematician, Floquet. This theorem states that for a given mode 
of propagation at a given steady-state frequency, the fields at one cross section differ 
from those one period away only by a complex constant. 

z, z,+d 
[-*- d 

aT 


Fig. 27. Definition of periodic planes (z\) and ( z\ + d) for the application of Floquet’s 
Theorem, and regions I and II for the formulation of the boundary value problem. 

The proof of the theorem lies in the fact that when a structure of infinite length is 
displaced along its axis by one period, it cannot be distinguished from its original self. 
(For a more formal proof, see, e.g., Landau and Lifshitz, Mechanics). 

Suppose we write 
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E(r, z, t) = F(r, z) e~^° z c** (4.47) 

where F(r, z) is a periodic function of z with period d, pod < it, ^-dependence is 
excluded, and propagation is assumed in the positive z-direction. Referring to planes 
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(z\) and (zi + d) in Fig. 27, expression (4.47) becomes: 

At z = z x E(r, z lt t) = F(r, z x ) e** , 

at z = zi+d E(r,zi+d,t) = F{r,zi + d)e~ j ^ Zl+d K^ t . 
But by definition 

F(r,zi + d) = F(r,zi) . 


(4.48) 


(4.49) 


Thus 


E(r, z\ + d,t) — E(r, z h t) e J/3 ° d , (4.50) 

quite independently of where z\ is taken, which indicates that the form of (4.47) indeed 
obeys Floquet’s Theorem. Since F(r, z) is periodic, we can expand it in a spatial Fourier 
series of the form: 

n==+o° 

F ( r > z ) = £ C a «( r ) e ~ J ^ 2 (4.51) 

n——oo 

where n is an integer which can take all positive and negative values. It can be shown 
that under rotational symmetry, the solution given in (4.47) becomes 


n=+oo 


£■(<•, *,()= Y. C a n Mk r „r) e'V-M 


n =—oo 


(4.52) 


where 


ik =00 + -j- , (4.53) 

, (4.54) 

C is the usual amplitude factor, and a n is the amplitude of the n th so-called space 
harmonic. 

_ There are several methods of determining the a n amplitudes of these space har¬ 
monics for a given structure. Before we mention them, we must make a number of 
observations regarding Eqs. (4.51) through (4.54). 
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The fact that the expansion contains all positive as well as negative integers 
arises from the fact that we have assumed an even structure (with reflection 
symmetry with respect to the origin, z = 0). 

The wave E z {r,z,t) which meets the boundary conditions of the periodically 
loaded structure is made up of space harmonic components which each have a 
different propagation constant (4.53). However, once /3q(oj) is known, all /? n ’s 
are known. 


Each space harmonic has a different phase velocity given by 

’** = • |4 ' 55) 

For any n < 0, the phase velocity of the space harmonic is negative, i.e., 
a particle would have to travel in the negative direction to remain in phase 
with it. The continuous lines in Fig. 26 all belong to a wave whose energy (or 
“group”) travels in the positive direction, and conversely for the dotted lines. At 
the u>o p intercepts, since the u — ft curves are simply translations of themselves 
by 2m?/d , all have the same slope corresponding to the same group velocity. 


In the definition of shunt impedance (3.8) and r/Q (3.14), the values of Vo/d 
that should be used are those of the space harmonic which is synchronized with 
the particle. In most cases, that is the fundamental, i.e., oq. 


From (4.52) and (4.54) it is seen that the terms a n Jo(k rn r) are not constant as 
a function of radius. If k^ n is negative, k rn becomes imaginary and Jq changes 
into Iq, which instead of falling off with r, increases. There is one notable 
exception which turns out to be extremely fortunate for the operation of linear 
accelerators. Indeed, when n = 0 and v p = c, /?q = k“ and k r o = 0, which 
makes Jq = 1 and the accelerating field independent of radius. As a result, 
acceleration and energy acquired are independent of the radial position of the 
synchronous particle. 


The actual variation of u> versus /?q depends on the relative diameters (2a and 
26) of regions I and II in Fig. 27. When 2a 26, the passbands are very 
narrow until, when 2a = 0, there is no more passband at all. Conversely, when 
2a —► 26, the passband becomes very wide until finally, we get back to the pure 
TMqi circular waveguide mode. 
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— The choice of the phase shift per period for the fundamental is generally made on 
the basis of where the highest shunt impedance per unit length can be obtained. 
Snapshots of electric field configurations for the cases where ftod = 0, tt/2, 
2x/Z and x are shown in Fig. 28. Somewhat loosely, these 0, tt/2, 27t/3 and x 
configurations are also called modes, not to be confused with TM, TE modes 
or superpositions thereof. Note that these field configurations have been chosen 
at such an instant in time that they can interchangeably represent a traveling 
wave or the maximum of a standing wave. As will be seen later, for traveling- 
wave accelerators of the disk-loaded type, the optimum choice is fad = 2xfZ 
(120°). For standing-wave accelerators, it is Pod = x or x/2, as shown in Sec. 
4.9). 



Fig. 28. Snapshots of electric field configurations for disk-loaded structures with 
various phase shifts per period (also loosely called “modes”). 


The u> — ft diagram can be determined exactly, point-by-point, by obtaining the 
resonant frequencies of stacks of n cavities, which have (n + 1) resonances, one for 
each value of ftod. For example, a stack of three cavities can resonate for ftod = 0, 
Pod = x/Z, ftod — 2x/Z and ftod — x. The frequencies can be obtained either 
experimentally or by means of computer programs such as SUPERFISH or LALA 
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which find the resonances and calculate the fields for standing waves. A conversion 
formula can then be used to obtain the fields for the traveling waves. The fields 
can also be found experimentally by means of bead perturbation techniques in cavity 
structures. An example of a field profile for the fod = 7r/2 “mode” is shown in Fig. 
29. 



Fig. 29. Example of approximate axial E^-field amplitude for /Sq d = n/2 “ mode.” 
This is a snapshot of a traveling wave at a given time. 

Yet another theoretical approach to obtain the fields is to set up general solutions 
in the two regions (I) and (II) shown in Fig. 27 (sometimes called “tube” and “slot” 
fields) and then to match these solutions along the common boundary. This approach 
was first used successfully in the early days (1948) by W. Walkinshaw, in England, 
long before the advent of the modern field-mapping computer codes. 

Once the fields are obtained, a simple Fourier-type analysis can be used to get the 
space harmonic amplitudes. The other structure parameters such as a, Q, v p , v g , r 
and r/Q can then also be calculated from the definitions and by the techniques intro¬ 
duced in Sec. 4.2. This entire subject is discussed in detail in a later chapter by E. 
Keil. 

Before we leave the subject of space harmonics and Brillouin diagrams, we should 
return briefly to the second passband shown in Figs. 23 and 26. As we hinted earlier, 
this passband has its origin in the second “normal mode” of oscillation of the cells mak¬ 
ing up the chain. From what we have learned about the lowest passband, the n = 0 
space harmonic of this higher band propagates in the negative direction ( doj/d/3 < 0). 
To call attention to this fact in Fig. 29, we have drawn this branch as a dashed line. 


60 



Because of this property, this passband is often referred to belonging to a backward 
wave, i.e., the fundamental has its phase velocity and group velocity going in opposite 
direction. 

The existence of this and higher passbands can be responsible for serious insta¬ 
bilities in the operation of linear accelerators, even though the RF sources driving 
them do not in principle supply power at these frequencies. What happens is that the 
short bunches of relativistic electrons can excite the structure at frequencies at which 
these higher modes are synchronous or quasi-synchronous with the bunches, whereby 
their amplitude can grow, either through feedback or cumulative interaction. This 
interaction can lead to deleterious effects such as beam deflection and ultimate beam 
breakup. Such effects can limit the achievable beam current in multi-bunch beams 
and cause emittance growth in single-bunch beams. These problems are treated in a 
later lecture on instabilities by A. Chao. 

4.8 Traveling Waves or Standing Waves? 

Until now, not much has been said about how the fields are established in the linac 
and where the rf power is fed into the structure. There are basically three ways in 
which this can be done. 

In the first, the structure is of the traveling-wave type, in which case the power 
is fed in at one end, propagates through the structure with some attenuation and the 
balance is absorbed in a load (Fig. 30a). Steady-state is reached when the structure 
is filled with energy after one pass. For a length £, the filling time is tp = 1/ v g . 

In the second, the structure is of the standing-wave type, in which case there is only 
one coupler, either at one end or in the middle of the structure (Fig. 30b). The fields 
build up through multiple reflections, the number of reflections of any consequence 
being determined by the wall losses. If the structure is long and the wall losses are 
high, the number of reflections that add any significant field is small, and vice versa. 
If the structure is matched to the power source, there is a so-called Qp or loaded Q 
associated with it which is equal to half the Qq associated with the wall losses. The 
maximum field build-up is attained with an exponential time constant tp — 2Qp/u 
or Qo/u, as defined in expression (3.16). In three filling times, the fields attain 95% of 
their ultimate value. One of the questions that must inevitably be asked is: aren’t the 
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reflected waves useless since they have their phase velocity in the wrong direction and 
hence do not deliver any net energy to the beam? The answer is that this is indeed true 
for all values of d tt. However, when fod = tt, the n = 0 space harmonic of the 
forward wave coalesces with the n = —1 space harmonic of the backward wave and 
both contribute energy to the beam if it is in synchronism. As a result, all standing- 
wave linacs are built in the so-called jr-mode, or in a variation thereof, as we will see 
in the next section. - - 
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Fig. 30. Linacs with (a) traveling-wave, (b) standing-wave and (c) traveling-wave with 
feedback structures and source power P e . 

The third method, shown in Fig. 30c, which was proposed in England in 1949 by 
R.-Shersby-Harvie and L. Mullett, is to feed the power left over at the output end of 
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the traveling-wave structure back to the input via a low-loss waveguide. At this point 
this power can be recombined at the correct phase with the source power by means 
of a bridge. In this scheme, the backward-wave problem is circumvented at the cost 
of considerable complexity. This scheme has few advantages over the standing-wave 
scheme and has been used rarely. 

The question of how one chooses between a traveling-wave structure and a standing- 
wave structure depends on the desired application and a number of other constraints 
such as the available rf power, energy and length, and the required pulse length and 
duty cycle. 

The following general comments and comparisons can be made: 

1. Since the standing-wave structure is generally built with a fixed coupling, it 
cannot be matched under all conditions. One can show that if the power source 
has a relatively short rf pulse length (measured in a few tp' s), at best 80% of 
the rf energy can be delivered to the structure during the pulse. The remaining 
20% is reflected to the source. This mismatch can under certain conditions 
be destabilizing to the source, causing frequency pulling, phase oscillations and 
fluctuations in output power. Furthermore, since the fields take several filling 
times (2Ql/lo) to build up to their final value, the beam will vary in energy 
throughout the rf pulse if it is injected after, say one filling time, and the total 
pulse is only a few filling times long. If the application calls for a precisely 
defined particle energy, this is then an unacceptable disadvantage. Therefore, 
in all such cases, the traveling-wave linac, which reaches steady-state in one 
filling time, is favored over its standing-wave counterpart. 

2. On the other hand, standing-wave structures such as the side-coupled or the 
disk-and-washer structure can be designed to have shunt impedances as much 
as 50% higher than traveling-wave structures at the same frequency. This ad¬ 
vantage also prevails for phase velocities in the range 0.3 < v/c < 0.8 where 
disk-loaded structures are relatively inefficient. If the pulse length is long or in¬ 
finite, i.e., the machine is CW, whatever energy is wasted during the filling time 
is negligible. Conversely, the traveling-wave structure, for a long pulse, always 
wastes in a load that power which is left at the output end. Thus, long-pulse 
proton linacs like the LAMPF accelerator, and all CW linacs use standing-wave 
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structures. Two examples of CW machines are the room temperature structures 
proposed for microtrons, and the superconducting structures which reach shunt 
impedances 10 4 to 10 5 higher than at room temperature. A third example is the 
structure used in all e ± storage rings. In this case, the standing-wave structure 
meets the extra requirement that it can accelerate particles of opposite charges 
in opposite directions. 

3. There is one case which falls between those described above. That is the case 
where one wants to obtain a high gradient in a short length with a source of 
relatively low power (e.g., a magnetron) and a pulse length of a few filling times. 
Typical applications of this case are the short medical accelerators and special 
sections for injectors and positron sources. The standing-wave structure here is 
at an advantage because the one-way wall losses are low and the large number of 
reflected waves can build up the fields to a high level. Were the reflected energy 
to be dissipated in an external load after a single pass, most of the power would 
be wasted. The only way to make the traveling-wave structure competitive in 
this situation would be to greatly increase the filling time of the short section 
by lowering its group velocity accordingly. This in turn would increase the 
energy density and the fields. However, to achieve this very low group velocity, 
one would have to reduce the iris diameter considerably. This would lead to 
beam transmission problems and also to difficulty with dimensional tolerances. 
Consequently, in these cases, designers are willing to live with the disadvantages 
of the standing-wave structure discussed earlier in order to achieve the required 
high energy within a very short length. 

4/9 Resonant Coupling 

There is one problem in the above discussion of standing-wave structures that still 
needs to be elaborated because it has led to important innovations. We have seen 
that using the jr-mode is indeed the way to make use of both forward and backward 
waves. The problem, however, is that if we resonate a structure of reasonable length, 
say 50 cavities, at u 0 p where 0 O pd = tt, it will contain 50 half-wavelengths at this 
frequency. Adjacent to it, however, there will be a frequency at which it will contain 
4fFhalf-wavelengths or for which fid = 49/50 n. Since the u — fi curve has zero 
slope at (3opX = d, the neighboring resonance at fid = 49/50 tt will be extremely 
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close to u> op , to the point where any small cavity mistuning, frequency modulation or 
beam loading will excite the neighboring mode (or modes). In the jargon, the ff-mode 
is said to be very unstable. If, on the other hand, u op were in the middle of the 
pass-band, i.e., at j3 op d — n/2, the slope of the oj — j3 curve would be greatest, and 
maximum frequency separation (or mode stability) would be achieved. This line of 
reasoning led E. Knapp at LANL to invent the side-coupled cavity structure in the 
early 1960’s. Figure 31 makes the invention understandable. Figure 31a shows the 
standing-wave structure resonating in the ;r/2-mode. Clearly, this configuration is 
inefficient compared to the 7r-mode since under steady-state conditions, every second 
cavity is unexcited and produces no acceleration. Figure 31b shows a modification 
where the excited cavities have been stretched and the unexcited cavities have been 
shrunk. Every cavity however, short or long, is still tuned to the same frequency. Thus 
the 7r/2 mode and its inherent stability are preserved. Finally, Fig. 31c illustrates the 
discontinuous step in the invention: the empty “useless” cavity has been placed on 
the side, off the beam line. From the rf point of view, the mode remains jt/ 2 but for 
the beam, it has become 7r-like. The side-cavities which do the coupling are empty 
but resonant, the feature which is responsible for the name of “resonant coupling.” 
Figure 31d shows the physical embodiment of the idea. The accelerating cavities are 
shaped and provided with nose cones for maximum shunt impedance, and the coupling 
cavities are staggered to reduce asymmetries introduced by the slots. 

Another way of looking at what has been done is to consider the oj — plots which 
accompany Fig. 31a, b and c. The double periodicity (in 31b) would normally lead to 
a break and a stop-band in the a > — /3 curve at the n/2 point: two branches are formed 
as in Fig. 26. However, through the process of resonant coupling, what is called “con¬ 
fluence” between the two branches is achieved, as shown in Fig. 31c, re-establishing 
the stability of the 7 t/ 2 mode. Thus in conclusion, at the cost of some extra complex¬ 
ity, the side-coupled cavity chain and a number of other structures based on the same 
principle have yielded the best of two worlds over a wide range of particle velocities: 
good efficiency of the TT-mode, and insensitivity to fabrication errors and beam loading 
of the ^/2-mode. The resonant posts or slugs used in the Alvarez structure shown 
earlier in Fig. 12 are another example of the same idea at lower frequency. 
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Fig. 31. Evolution of n /2-mode iris-loaded structure to jr-mode side-coupled structure. 
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5. Electron Linear Accelerators 


In this last section we will cover in very broad terms a few of the major design 
and operating features specific to electron linacs. A number of references will be made 
to the SLAC linear accelerator as an example. For more details, the reader should 
consult the book on the Stanford Two-Mile Accelerator, R. B. Neal, General Editor, 
W. A. Benjamin Inc., 1968, hereafter abbreviated S.T.MA- . asjwell as L.A. referred 
to earlier. 

5.1 The Main Elements 

The main elements of an electron linac are shown in Fig. 32. The electron gun is 
a triode consisting of a thermionic cathode, a wire-mesh control grid and an anode. 
The potential difference between the cathode and the anode is typically 80 kV which 
gives the electrons a velocity of c/2. The grid is normally held at a voltage somewhat 
negative (~ 50 V) with respect to the cathode so that the gun is biased off except 
when a beam pulse is needed. To trigger the gun, a positive pulse, typically of a 
few hundred volts, is applied for the desired length of time. Fast transistor pulse 
circuitry has improved in the last few years to the point where beam pulses as short as 
one or two nanoseconds are now achievable. Most electron linacs operate with beam 
pulses in the 1-10 nsec range at repetition rates up to 500 and even 1000 pps. The 
continuous stream of electrons out of the gun first passes through a prebuncher. This 
is a single low-Q re-entrant cavity supplied with typically ~ 2 kW of rf power at the 
fundamental frequency. The function of this prebuncher is to velocity-modulate the 
electrons with peak gap voltages up to 10-15 kV so that early particles are slowed down 
and late particles are speeded up. As a result of this process, longitudinal compression 
or bunching results. 

At a distance of 30 cm downstream, one can get on the order of 70% of the electrons 
into a phase interval of 70 electrical degrees. At this point, a refinement which is not 
shown in the figure is sometimes used to discard a large fraction of the bunches. This 
is a transverse electric field rf beam chopper which can act as an additional gate (say 
for single bunches every 12.5 nsec if the chopper frequency is 40 MHz) by deflecting all 
bunches except the ones that go through at zero-crossing. Yet another variation is to 
combine the prebuncher and the chopper into one rf cavity at a subharmonic frequency 
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of the accelerator. This scheme then produces a train of more widely spaced and longer 
bunches with a potential of packing more charge into each of them. Combined with 
a fast grid pulser, it can generate a single bunch. This scheme is being used for the 
SLAC Linear Collider (SLC) injector. 
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Fig. 32. An electron linac (or the beginning of one!) 


The next element in Fig. 32 is the buncher. At SLAC this is a 4-cell-long, 27r/3- 
mode, traveling-wave section with a constant phase velocity of 0.75 c and a maximum 
power input of 2 MW . The bunches emerge from it with a kinetic energy of about 250 
keV in a phase interval of 30°. Bunching takes phase here because on the average the 
entering bunches are made to ride ahead of the crest and early electrons at first get 
slowed down while later ones get accelerated. Being slower than the wave when they 
enter the section, they then all rise towards the crest and gain energy. In addition 
to this average motion, they perform approximately one quarter-wave oscillation with 
respect to the central electron: the buncher acts like a quarter-wave impedance match¬ 
ing transformer in the e — <I> space of Fig. 6, between the prebuncher on the one side 
and the accelerator section on the other. The accelerator section is a traveling-wave 
structure with phase velocity c. In the SLAC case, as will be described later, it is an 
S-Band, constant-gradient, 27r/3-mode, 3-meter long structure. The 30° long bunches 
are caused to enter it around the field null ahead of crest, and they asymptotically 
approach the crest of the wave as they are accelerated. Unlike the protons at lower 
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velocities, they do this monotonically without oscillations. A simple expression for this 
asymptotic motion will be given in sub-section 5.3. At the output of the first accelera¬ 
tor section, the electron bunches are 5° long and their energy is typically on the order 
of 30 MeV for an rf power input P of 16 MW. This power is typically supplied by a 
high power klystron which in turn is pulsed by a high voltage modulator. RF power 
is also supplied to the buncher and to the prebuncher. 

In addition to these accelerating and bunching fields, the electrons are acted upon 
by their own space charge forces, radial rf fields from the structures, stray magnetic 
fields, and in case of very high currents, Wakefields which they themselves induce and 
leave behind. For all these reasons, steering, degaussing and focusing devices, not 
shown in Fig. 32, are also incorporated into the design of the linear accelerator. As 
an example, in the SLAC injector, there are two thin magnetic lenses to focus the 
beam. The first lens is located downstream of the gun, and the prebuncher gap is 
placed at the beam waist formed by this first lens. The second lens is located after the 
prebuncher and produces a second waist at the entrance of the buncher. These waists, 
i.e., points of minimum beam diameter, are created at those points where one wants 
to minimize the effects of radial fields which generally grow linearly with radius. 

Along the buncher and the accelerator section there is almost always a solenoidal 
DC magnet field of a few kilogauss. The required focusing action of the axial magnetic 
field is obtained because there are radial DC fields where the beam enters the solenoid. 
This produces an azimuthal force causing the beam to begin to spiral. The azimuthal 
velocity and the axial field produce the desired inward force. 

The linear accelerator which has been described so far either ends at this point, in 
which case the beam is extracted for whatever use it is destined; or it just constitutes 
the injector or beginning of a longer, multisection machine such as the Orsay, Saclay, 
Kharkov, KEK, MIT, NIKHEF, SLAC and many other linacs in the world. The design 
that has been described is fairly standard. The only major variance is in the buncher 
which in a number of short accelerators is designed with a tapered phase velocity which 
gently forms and “escorts” the bunches up to the velocity of light. The buncher and 
the accelerator in this case are constructed in one section powered at the upstream 
end. In the low energy (4-6 MeV) medical linacs, the rf source is a 2MW magnetron 
and the structure, as discussed in Sec. 4.8, is generally of the standing-wave type, 
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under one meter of length. Solenoidal focusing can sometimes be avoided by using the 
alternating-gradient focusing effect of the TT-mode radial fields in each cavity. 

Referring back to Fig. 1, it should be noted that in all these machines, when 
the electrons have acquired a kinetic energy of 3 MeV, they have already reached a 
velocity of 0.99 c, i.e., they are totally relativistic. Hence, the beam dynamics “drama” 
described earlier for protons in Secs. 3.5 and 3.6, which takes place over hundreds of 
meters in a machine such as LAMPF, is played in an electron linac in the first meter! 

The longer, higher energy electron linacs simply use many separately powered 
sections. The trade-off between rf power and length is generally resolved on the basis 
of economics, technological factors, the availability of land and the cost of AC power 
to operate the machine over a number of years. In addition to the rf structures, 
klystrons and modulators, these longer machines use periodic focusing systems in the 
form of FODO arrays, quadrupole doublets or triplets. Of course they also have their 
associated vacuum, mechanical support, cooling, instrumentation and control, and 
other support systems. For positron operation, they generally have a tungsten target 
which is inserted at a point where the electrons have achieved a few hundred MeV, 
or at SLAC 5-10 GeV. The positrons are generated in a wide forward cone which 
must be focused by a high field solenoid so that the beam can be contained within the 
transverse admittance of the subsequent accelerator sections. 

All these subjects are described in much greater detail in L.A. and S.T.M.A. 

5.2 Structure Design, Choice of Parameters and Energy Gain 

Having described the main elements of a “generic” electron linac, we will now re¬ 
turn once more to the structure and to some of the practical choices involved in its 
design. The first choice is to select the general type of structure. In what follows we 
will concentrate on the simple traveling-wave disk-loaded structure. Only a few com¬ 
ments will be made on standing-wave structures for the sake of comparison. The most 
critical parameters to be chosen are the frequency, the phase-shift per cell, and the 
attenuation per section. On the attenuation, there is in addition a sub-choice between 
two possibilities. Either the attenuation is kept constant as a function of length in the 
section and the structure is uniform, in which case the power decays exponentially: the 
structure is called “constant impedance.” Or the attenuation is increased as a func- 


70 


tion of length by tapering the dimensions of the structure so that the power decays 
linearly and the field is kept constant: the structure is called, “ constant gradient .” To 
understand the effect of these parameters, it is best to calculate the electron energy 
gain in terms of the shunt impedance per unit length, redefined once more as 


_ E 2 
r ~ -dP/dz 

where E is now the electric field amplitude of the fundamental, 
harmonic. 


(5.1) 


synchronous space 


The other quantities of interest defined earlier are: 
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(5.5) 

from which one can get the relation 
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For a section of length £, the kinetic energy gain per section for a synchronous electron 
riding on the wave crest is then 

AE K = e f * E{z) dz . 

(5.7) 

For the constant-impedance case, a is not a function of 

z and from (5.1) we get 

E 2 (z) = 2ar P{z) 

(5.8) 


whereby the integration of (5.7) leads to 

A E k = e ( 1 —(2r) 1/2 (P 0 ^) 1/2 (5.9) 
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where Fq is the input power to the section and 


t — at 


(5.10) 


For the constant-gradient case, the form of a(z) that leads to a constant E(z) must 
be determined. To do this in closed form, we will assume that we can neglect the 
variations of r and Q that occur if the structure dimensions ace modified to adjust 
0 ( 2 ), and that, in fact, the only quantity that changes rapidly is the group velocity in 
(5.6). Since from (5.1) 

E 2 (z) = -/-^ , (5.11) 

if one wants to make E(z) constant, dP(z)/dz must be kept constant. Then redefining 
(5.10) as 



the linear decaying profile of P(z) must be given by 


P( 2 ) = P„ 1 - f(l - e~ 2r ) 


whereby 


and 


<*(z) = 


=-27- 
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Vg(z) = 
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(5.12) 


(5.13) 


(5.14) 


(5.15) 


We can see that to get a constant-gradient profile, we must let the group velocity 
decrease linearly with the same slope as the power. This is obvious if we refer to (5.5) 
since then w s t is also constant. We can now obtain the kinetic energy gain for the 
constant-gradient section. Expression (5.7) is simply A Er = eEt and since, from 
(5.13) 

f =-^(l-e- 2 ') , (5.16) 
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using (5.11) we obtain 


AE K =e(l-e~ 2T ) l / 2 (P 0 re) 1 / 2 . (5.17) 

Using expressions (5.9) and (5.17), we can now draw some very general conclusions. 

Both expressions depend on the factor (Pq r t ) 1 / 2 . The trade-off between Pq and t 
mentioned earlier is immediately apparent. For maximum energy gain, one wants to 
maximize r and the factors of r. From Fig. 33, we see the behavior of r as a function 
of the number n of disks per wavelength. The solid curves are calculated for an array 



n 


Fig. 33. Shunt impedance per unit length versus number n of disks per wavelength 
for various thicknesses t. 

of pillboxes of diameter 26 with zero-iris diameter (2a = 0). For this reason, they are 
somewhat higher than the experimental values shown by the dashed curves. Both sets 
are taken for three values of disk-thickness t, v p /c = 1 and / = 2856 MHz. It is 
seen that a broad maximum is obtained for n = 3 disks/Xo, i.e., for the 2jr/3-mode 
discussed earlier. For this reason, this mode was selected for the SLAC linac. The 
thicker value of t (0.230 n ) was selected for greater mechanical strength and because 
it permits a greater radius of curvature on the disk edge, which decreases the risk of 
electrical breakdown. 
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The factors of r in Eqs. (5.9) and (5.17) are compared in Fig. 34. It is seen that for 
small t, they are about equal and that for r ~ 1, they are close to 0.9. The value of 
t = 0.57 was selected for the SLAC linac as a broad compromise between maximizing 
energy and minimizing the effect of beam loading (see next sub-section) and the filling 
time (tp = 2 Qt/cj). The constant-gradient design was selected because of a number of 
advantages such as lower ratio of peak-to-average field and uniform power dissipation 
per unit length. It also turned out ex-post-faeto that tbe constant-gradient structure, 
due to its non-modular construction, is less susceptible to cumulative beam breakup 
because the effective interaction length with the synchronous transverse deflecting 
mode is foreshortened. 



T 

Fig. 34. Normalized energy gain versus r. 


; The choice of rf frequency is also made on the basis of a broad set of criteria. Table 
4, which is taken from S.T.M.A. , shows the frequency dependence of some of the key 
parameters. The single most important quantity is the shunt impedance r which varies 
as Z 1 / 2 . Since Q falls as / -1 / 2 , ru/Q which is a measure of E 2 /w s t varies as / 2 . This 
is obviously correct since the volume per unit length of structure varies as /~ 2 , all 
transverse dimensions varying as / -1 . Thus, any linac design for which one wants 
to maximize the beam energy and in which the stored energy is discarded after each 
pulse, will be favored by high frequency. On the other hand, the frequency cannot be 
increased indefinitely because the decreased iris diameter (2a) eventually becomes a 
limitation, both from the point of view of beam interception and transverse Wakefields 
which scale roughly as (2a) -3 . Also, at lower frequency, the larger energy stored makes 



the structure less sensitive to beam loading. Finally, there are factors which, although 
not as fundamental, can be just as important such as the availability of high power 
rf sources as a function of frequency, and the absolute machining tolerances. For all 
these reasons the SLAC frequency was set at 2856 MHz or 10.5 cm wavelength. Most 
other linacs in the world are built at the same frequency or nearby, in the vicinity of 
3000 MHz. The final dimensions of the SLAC sections are shown in Fig. 35, together 
with the linear profile of the normalized group velocity^ Vgjc) and the resulting shunt 
impedance. As we see, r(z) has an average value of 57 Mf)/m and it varies slowly with 
z, which makes our original assumptions acceptable. Figure 36 is a photograph of the 
final 3.05 meter-long SLAC constant-gradient section with its input and output 
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Fig. 35. Variations of 2b, 2a, v g fc and the shunt impedance r (corrected for the 
fundamental space harmonic) as a function of cavity number along SLAC 3.05 meter 
constant-gradient section for t — 0.230 in. 
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Table 4. 

Frequency Dependence of Principal Linear Accelerator Parameters 



Frequency 

dependence 

Frequency 

preference 


Parameter 

Higffi 

Low 

Notes 

Shunt impedance per unit length (r) 

pn 

X 


a 

RF loss factor ( Q) 

r x ' 2 


X 

a 

Filling time ( t r ) 

r 3n 

X 


a, b 

Total RF peak power 

r ui 

X 


a, b, c 

RF feed interval (/) 

r 3 ' 2 


X 

a, b 

No. of RF feeds 

pi 2 


X 

a, b, d 

RF peak power per feed 

r 2 

X 


a, b, c 

RF energy stored in accelerator 

r 2 

X 


a, b, c 

Beam loading ( —dV/di) 

pn 


X 

a, b, d 

Peak beam current at maximum 
conversion efficiency 

r xn 


X 

a, b, c, f 

Diameter of beam aperture 



X 

a 

Maximum RF power available from 
single source 

r 2 


X 

e 

Maximum permissible electric field 
strength 

p,2 

X 


9 

Relative frequency and dimen¬ 
sional tolerances 

P' 2 

X 


a, b 

Absolute wavelength and dimen¬ 
sional tolerances 

r U2 


X 

a. b 

Power dissipation capability of 
accelerator structure 

/- 


X 

a, b, d 


Notes: 

a. For direct scaling of modular dimensions of accelerator structure. 

b. For same RF attenuation in accelerator section between feeds. 

c. For fixed electron energy and total length. 

d. For fixed total length. 

e. When limited by cathode emission. 

f. When limited by beam loading. 

g. Approximate; empirical. 
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Fig. 36. SLAC 3.05 meter-long constant-gradient accelerator section. 

couplers. Using Eq. (5.17), the reader should verify that the energy obtained per 
section, taking into account ~ 10% feed losses, is approximately 

AEk = 10 e \JP\fw MeV . (5.18) 

If the power of a 36 MW klystron is split four ways, as it is at SLAC, then 

A Ex = 30 MeV/section 

or 120 MeV/klystron. 

In comparison, a matched standing-wave structure yields a steady-state energy of 

A E k = e (iV) 1/2 (5.19) 

with values of r that can be 50% higher. Thus, in principle, energies that are 25% 
higher could be obtained for the same values of Pq and l. In actual fact, because of all 
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the difficulties listed under Secs. 4.8 and 4.9, such standing-wave structures are only 
built for the special applications discussed earlier. 

5.3 Beam Current, Emittance, Bunch Length and Beam Loading 

In the preceding sub-section, not much has been said about the beam itself, its 
characteristics, its behavior and its effect on the structure. The subject of beam 
dynamics in an electron linac is again so vast that we cannot hope to do it justice in 
the space allotted here. Thus only a few aspects of the problem will be touched upon. 

A bunch of electrons can generally be characterized by the total number of particles 
contained in it, the transverse emittance, the bunch length and the energy spectrum. 
In addition, the beam, which is generally made of a train of bunches, has a pulse 
length and a repetition rate (unless it is CW), and the notion of beam power can be 
associated with it. We will now discuss these various properties, illustrating them with 
some practical examples. 

Total number of particles per bunch. Typically, in an electron linac, the total num¬ 
ber of electrons per bunch is in the range 10 ^ - 10 11 . For example, the SLAC linac, 
when it accelerates a beam of 50 mA peak current within a 1.6 //sec pulse, produces a 
train of ~ 5000 bunches and a total of 5 X 10 11 electrons per pulse. Thus, the number 
of electrons per bunch in this case is 10 8 and the beam energy per pulse at 25 GeV is 
2000 joules, or 720 kilowatts at 360 pulses/sec, i.e., close to one megawatt! The goal 
of the SLC project at SLAC is to obtain single bunches of electrons and positrons of 
5 X 10 10 particles, i.e., of 500 times greater intensity. Such single bunches with ~ 8 
nanocoulombs of charge at 50 GeV will carry an energy of 400 joules, or 72 kilowatts 
at 180 pulses/sec. 

By contrast, a medical linac at 5 MeV, 200 mA peak current and a duty cycle of 
10 -8 can deliver to its x-ray target a power of one kilowatt. 

Transverse emittance. Assuming a bunch of Gaussian or other not too unreasonable 
shape, one can define transverse emittances a x i and where a x and are 

the rms transverse dimensions of the bunch and <^ x and a y are the corresponding 
rms angular dimensions; for a single particle J — p x /pz = Px/7/frnoC and = 
pjjp z — pz/'i pniQC. The p x and py s refer to transverse momenta and 7 is the usual 
(1 —/3 2 ) -1 / 2 . Often, the above quantities are defined as /? 7 < 7 Z < 7 *. and P^VyC^y, or simply 
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7 (Tx<^x aQ d 'IVyo’y w h e n fl —*■ 1- These quantities have the advantage of being invariant 
with energy. They are measured in many different units of which some typical ones 
are mgc • cm or (MeV/c) cm. Since the second unit is twice as big as the first, the 
emittance measured in (MeV/c) cm is half that in moc • cm. In the x, p x and y, p y 
phase-space planes, Gaussian beam envelopes are represented by ellipses and one often 
quotes the above quantities multiplied by jt to denote emittance areas (see Fig. 37). 
According to Liouville’s Theorem, these areas are invariant provided that only non- 
dlssipative forces act on the particles. 


Fig. 37. (a) Logarithmic orbit due to 
initial deflection 9q at zq showing 
contracted length L; (b) Rotation of 
transverse phase-space ellipse from 
injector to accelerator output. 
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As an example, a recent transverse emittance measurement for a single bunch of 
5 X 10 10 electrons out of the proposed 40 MeV SLC injector gave ~ 7.5 X 10 -3 tt 
(M eV/c) cm. This means that if this beam is focused to a waist of 0.75 cm radius, its 
angular divergence is 0.25 mrad. If there were no transverse emittance growth along 
the accelerator, at say 40 GeV, the angular divergence would shrink to 0.25 jirad for 
the same 0.75 cm waist (see proof below). Even though the above emittance from 
the injector is quite small, the SLC requires that the emittance at the final focus be 
no greater than 1.5 X 10 -3 ir (MeV/c) cm or five times smaller. This explains why 
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damping rings to “cool” the beams downstream of the injector had to be added in the 
SLC design. 

The transverse beam dynamics in the injector of a linear accelerator are compli¬ 
cated by the presence of space charge forces and radial defocusing forces due to waves 
and particles which are not yet moving at the velocity of light. We have seen earlier 
from Eq. (3.43) that as (3 S —► 1, the magnetic and electric Jorces-from a synchronous 
wave cancel each other. Thus, leaving out all other spurious forces, the transverse 
equation of motion of an electron in the x- (or y)-plane past the injector is given by 


dp x 

dt 


= 0 


(5.20) 


but since dz = c dt, we can rewrite this equation as 


Integrating, we obtain 


d ( dx\ 
dzV dz) ~ ' 


dx 

7 — = const 
dz 


(5.21) 


(5.22) 


which confirms the correctness of the assumption we made above. 

Since under these conditions the quantity 7 dxfdz is conserved along the length of 
the linac, we can make it equal to 70 0o at zq, say at the output of the injector. Then 
a second integration leads to 

/ x r z dz 

dx = 7o0o / — (5.23) 

X 0 J Zq 7 

and if we assume that the energy grows linearly as 7 = 7*2 where 7 * is constant, then 


x — xq = ~ 0q in— (5.24) 

7 zq 

or 

z 

x — xq = 0q zq in — . (5.25) 

_ Zq 

Looking at Fig. 37a, we see that an electron starting at zq with energy 70 (where 
actually zq is defined as 70 / 7 *), x — xq and angle 0 q, will end up at 2 with transverse 
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displacement zo + 0o 2 O^ n z f zq. It is as if the length z — zq has been contracted to 
L = zq£h z/zq. Thus if zq = 3 m, z — 3000 m, q/ = 26.66/m, the 2997 m length 
looks like it has been contracted to 20.7 m. The orbits instead of being linear are 
logarithmic. 

It is instructive to follow the evolution of the phase-space ellipse from its waist 
at the output of the injector (40 MeV) to its final configuration at the output of the 
accelerator (40 GeV). Since p x cannot grow beyond (po)max and the two (xQ) max points 
represent rays which are parallel to the z-axis and therefore cannot be deflected by 
an increase in energy, the ellipse has no choice but to get elongated between the two 
(po)max horizontal lines while remaining fixed at the two (zo)max points. It can be 
shown that as the ellipse gets elongated, the (x) max coordinate is given by 

X max — (^olmax "b ^0^0 ~J (5.26) 

where 0q = ( Po)max/PzO • In the example we chose, Xq = 0.75 cm, 9q — 0.25 mrad, 
zq(ti z/zq — 20.7 m, we find that x max = 0.91 cm. We see that the radius of the beam, 
in the absence of any transverse emittance growth, would barely need any focusing to 
stay within the aperture of the accelerator irises. In actual fact, however, emittance 
growth does take place because of the effects discussed earlier, and a FODO array or 
other quadrupole focusing system is indeed required. The description of these systems 
is beyond the scope of these lectures. 

Bunch length. In sub-section 5.1, we have described qualitatively how the bunches 
are forfned in the injector of a linear accelerator. Following the electrons from the gun 
through a prebuncher and a buncher is a complicated problem which is best done by 
a computer. However, when the particles finally enter the v p — c constant-gradient 
accelerator section, there is a fairly simple expression for the asymptotic phase <^>oo of 
an electron of velocity t>o = A) c entering with phase <f>o relative to the wave: 


Unlike for protons, the phase origin is taken at the field null, 90° ahead of crest 
(Fig. 38) and a = e FAq/woc 2 . * When the factor 2jr/a f(l — A))/(l + A))] 1 / 2 is 
made equal to 1 through the proper choice of a (a = 3.63 if A) = 1/2, or a = 2.37 
if /3 0 = 3/4), then an electron, entering the accelerator at 4>q — 0, asymptotically 
approaches the crest, <j ><*> = —tt/ 2. For small phase extents ±^o around <f >o = 0, it 
can be shown that 




i.e., all other electrons end up behind the crest. For example, if —1/4 < <f >q < 1/4 
radian, then — [tt/ 2 + 1/32] < <^oo < —tt/ 2 radian, i.e., a 30° bunch ends up in less 
than a 2° bunch. Typically, at SLAC, bunches over a wide range of charge have 
been measured to contain most of their charge within 5°. The asymptotic charge 
distribution is not quite Gaussian, even if that of the entering bunch is (which is 
generally not the case either). For the SLC, it is planned to have most of the 5 X 10 10 
electrons within 20°. 


Fig. 38. Asymptotic bunching process in 
v p — c constant-gradient accelerator 
section with value of a optimized 
for entrance conditions. 



Energy spectrum due to bunch length and beam loading. The last beam property of 
great interest to linac builders and users which we will consider is the energy spectrum, 
namely, how monochromatic is the beam? Typically, one wants to have all the electrons 
within an energy spectrum width of 1% or less. There are two effects inherent to how 
an electron linac works which bear on the spectrum. The first one is the bunch length 
which by its very nature implies that all electrons within it do not receive the same 
energy, either because they are not on the same part of the wave or because of their 

*Not to be confused with the attenuation/unit length a defined earlier! 
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own wakefields. The second one is beam loading from bunch to bunch due to energy 
removal from the structure by each bunch: this effect can reach steady-state after a 
certain number of bunches. In addition to these two inherent effects, there are a number 
of other technological effects which can be detrimental to the spectrum width such as 
poor injection conditions, phase and frequency modulation in the master oscillator, 
poor amplitude and phase stability of the rf sources, poor accelerator temperature 
control and many others. We will not examine these. 

In a multi-section linac with individually phase controllable klystrons, it is possible 
to phase each tube so that a bunch of phase length A0 is centered on the crest. Then 
it is easy to show that 


6 A E k , A <f> (A <t>) 2 

7 |7 - 1 COS —— --- 

A E k 2 8 


(5.28) 


or ~ 0.12% for A0 = 1/10 radian. The reader can convince himself that this is the 
best that can be done for a low current bunch in which wakefields do not play a role. 
Indeed, this effect results from the shape of the cosine wave and nothing can be done 
about it, not even rocking the bunch from one side of the crest to the other, which 
is sometimes suggested! This would only work if the top of the wave were a triangle 
instead of a curve. On the other hand, if the bunch is highly charged as in the SLC, the 
energy decrease due to the wakefield effect of the head on the tail can be compensated 
for by locating the bunch ahead of crest at the proper phase, which to first order can 
cancel the decrease if the slopes can be matched. 


- Coming back to the low current case, it is often asked what the effect of misphasing 
individual klystrons is on the beam. The answer is that to first order, misphasing 
individual klystrons only decreases the maximum reachable energy but does not affect 
the spectrum. Indeed, it can easily be shown that the (A<£) 2 /8 rule still applies as long 
as on average, the bunch travels centered on crest. 


The beam loading effect, which has been referred to several times earlier, comes 
about because the rf energy stored in the linac structure is reduced as successive 
bunches get accelerated and extract energy from it. The effect goes through a tran¬ 
sient phase until a new lower equilibrium energy level is attained. In a traveling-wave 
structure, it can be seen intuitively that each bunch that traverses it sets up instanta¬ 
neous rf field “packets” in each cavity at the velocity of light (except in the buncher). 
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These packets then travel down the structure in phase at the group velocity. If the 
train of bunches continues, there is a build-up of fields versus z, even though the pack¬ 
ets get attenuated with distance. Steady-state is reached when the packet generated 
in the first cavity reaches the last. As it happens, this is exactly the definition of the 
filling time tjr. Thus if we inject a train of bunches after one filling time, i.e., when 
the section if filled, it takes an additional filling time before steady-state is reached. 

We will give a short derivation of the steady-state case and let the transient case 
be understood intuitively. If a peak current of intensity i is assumed, our power loss 
equation (5.4) must be written as 


dP 

dz 



{Ei)beam i 


(5.29) 


assuming that the electrons are synchronous and ride on the crest. Then, differentiat 
ing the expression 


E 2 — 2a rP 


(5.30) 


with respect to z, we get 


dE da dP 
E -j— — rP — + ar — 
dz dz dz 


(5.31) 


Substituting (5.29) for dP/dz in (5.31), —E 2 /r for ( dP/dz) wa u and £’ 2 /2ar for P, we 
obtain the general differential equation 

dE f 1 da\ 

Tz = -“K 1 - Tz) -"" ■ ,5 ' 32) 

It follows that for the constant-impedance case 


dE 

dz 


—a E — art 


(5.33) 


and for the constant-gradient case, replacing dafdz by its value given by (5.14), 

dE , 

rf7 = ~ a " ■ ,5 - 34) 


The solutions are: 
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Constant-Impedance 


E = E 0 e~ az - rt(l - e~ az ) (5.35) 

where Eq = (2a r Pq)^! 2 , 

Constant-Gradient 

- - E = Eo + t ^tn\l-~ (l-e“ 2T )] (5.36) 

where 

B 0 = (2a 0 rP 0 ) 1/2 = [~(l-^ 2r )] !/2 • (5.37) 

The energy equations are then obtained by simple integration, as in (5.7): 
Constant-Imp edan ce 

A E K = <P 0 ir) 1 ! 2 (2 r) 1 / 2 —~ - eiri^L - -1 (5.38) 

Constant-Gradient 

A^=e(P 0 <r) 1 / 2 (l-c- 2 ') I/2 -^(l- i ^^ ? ) • (5.39) 


We can see that when i = 0, expressions (5.38) and (5.39) collapse back to (5.9) and 
(5.17) respectively. The reader can verify that for the SLAC conditions used for (5.18) 
and t = 50 mA, r = 0.57, the energy reduction per section due to beam loading 
under steady-state conditions is 2 MeV out of 30 MeV or 6.7%. The steady-state field 
profiles along a section given by (5.35) and (5.36) are sketched in Fig. 39a. The shaded 
areas show the reduction due to beam loading. The energies are the integrals under 
the respective curves. Figure 39b gives the evolution of energy as a function of time. 
The interval between 0 and tp shows the energy increase as the section is filling, the 
second interval between tp and 2 tp shows the transient decrease due to injection of 
the beam at t — tp , and the third interval (f > tp) shows the steady-state energy. 

The beam loading effect, as can be seen, is inherent to the operation of the linac. 
Once steady-state is reached, it has no effect on the spectrum width, unless the current 
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itself fluctuates. If the beam pulse length itself is only on the order of 2 tp, which is 
the case at SLAC, half the beam pulse is in the transient state. The remedy which 
is often used to reduce this overall spectrum broadening effect is to inject the beam 
somewhat earlier than at time tp so that the first bunches see their energy reduced 
as well. With a large enough number of klystrons, it is possible to stagger the onset 
of the rf in successive sections so that the two effects cancel each other down to less 
than 1% (Fig. 39c). Thus, the steady-state reduction cannot be'avoided but at least 
ifs effect on the spectrum can be minimized. 
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Fig. 39. Effect of beam loading: (a) on £-field profiles in sections, (b) on energy A 
delivered versus time, (c) on energy A Ep; if beam is injected before t — tp. 
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In standing-wave accelerators, beam loading is of course also present. However, 
since the beam-induced waves are caused to bounce back and forth just like the ex¬ 
ternally supplied waves, the reduced field profile takes several filling times to reach 
steady-state. 

5.4 Recent Developments and Future Challenges 

During the past few years, progress in the design t>f linear "electron accelerators 
has been made on a broad front. Industrial manufacturers have specialized in the 
production of a large number of medical and radiographic linacs. Their contribution 
has been strong in the areas high rf-to-beam energy conversion efficiency, high shunt 
impedance standing-wave structures, minimum external focusing, good packaging, and 
in some specialized machines excellent emittance characteristics through the design of 
very high voltage electron guns (up to ~ 400 kV). In some cases, linacs have been 
built to produce beam pulses in the 10-20 nsec range with currents up to 20 A. In 
the'universities, the accent at laboratories such as MIT, ALS (Saclay) and NIKHEF 
(Amsterdam) has been on linacs with duty cycles up to 10%. The challenge there 
has been to obtain reasonable currents (~ 20 mA) for long pulses (10-40 //sec) with 
excellent emittances and without the onset of beam breakup. 

At laboratories in the U.S. such as HEPL and SLAC (Stanford), University of 
Illinois and Cornell, in Germany at Karlsruhe and Wuppertal, in Switzerland at CERN, 
and in Japan at KEK, considerable effort has gone into investigating superconducting 
structures for linacs, microtrons and storage rings. The hope to build CW linacs 
with accelerating gradients up to the theoretical limit of 30 MV/m in niobium has 
not been realized so far. On the other hand, much progress has been made towards 
understanding and overcoming some of the major obstacles in the field such as the onset 
of multipactoring and thermo-magnetic breakdown due to impurities in the niobium. 
Thus, with a few tens of watts of power, several laboratories have been able to obtain 
gradients of 2-4 MV/m fairly routinely under CW conditions. The entire subject of rf 
superconductivity in its application to accelerators is dealt by M. Tigner, elsewhere in 
this book. 

_Other institutions have built or used a variety of dedicated low duty-cycle linacs 
for injection of electrons (and positrons) into storage rings. This category includes 
several laboratories such as Cornell, Frascati, Orsay, DESY, Novosibirsk, the Photon 
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Factory at KEK, and in the future the LEP injector at CERN and two injector linacs 
being planned in the P.R.C. (Beijing and Hefei). A large 4 GeV linac, with 200 mA 
peak current and high repetition rate (1000 pps) is presently being designed as an 
injector to a stretcher ring for the National Electron Accelerator Laboratory (NEAL) 
at Newport News (Virginia). 


Another area where innovation has taken place is in the field of rf energy compres¬ 
sion. The attempt here has been to gain electron energy without increasing average 
or even klystron peak power. The basic scheme, invented at SLAC by P. Wilson, D. 
Farkas and H. Hogg, has been called SLED, an acronym for SLAC Energy Develop¬ 
ment. (For a general reference, see for example, Z. D. Farkas, H. A. Hogg, G. A. Loew 
and P. E. Wilson, Proceedings of the DC^ International Conference on High Energy 
Accelerators, SLAC (1974), page 576). The principle is illustrated in Figs. 40 and 41. 



Fig. 40. The SLED principle. 


Referring to the top of Fig. 40, we see that under normal non-SLED operation, the rf 
drive pulse which is amplified by the klystron is directly transmitted to the linac. The 
SLED system, which is shown at the bottom of Fig. 40, has two major components: 
aT80° fast phase shifter on the drive side of the klystron and two high-Q (Qo — 
100,000) cavities on the output side of the tube with a 3 db coupler connected as 
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shown. During the first part of the pulse, the phase of the rf drive signal is reversed 
and the rf cavities fill up with energy with that phase. Because of the 3 db coupler, no 
energy is reflected to the klystron. The fields emitted by the cavities add algebraically 
with the fields reflected by the cavity coupling irises, and the power flows toward the 
accelerator. Exactly one linac filling time (~ 0.8 //sec) before the end of the pulse, the 
phase of the drive signal is flipped back by 180°: as a result both cavities discharge 
their energy into the accelerator. The discharging pulse amplitude is increased by the 
fact that the klystron pulse adds to it. The output pulse length is tailored exactly to 
fill the accelerator sections at one instant of time. 

SLED OPTIONS 



I 

n 

PULSE LENGTH 

2.5/as 

5/as 

EFFECTIVE POWER 
GAIN 

2 

3.15 

ENERGY GAIN 

1.4 

1.78 

REPETITION RATE 

360 pps 

180 pps 



Fig. 41. Comparison of two rf pulse length options for SLED showing how a longer 
charging time increases the output pulse amplitude. 

It can be seen intuitively that the longer the time to charge the cavities, the more 
energy there is available at the time of the discharge. On the other hand, there is a 
limit to this process because of the copper losses in the cavities. Figure 41 shows two 
rf pulse length options in what is obviously a continuum. Much beyond 5 //sec, one 
reaches a point of diminishing returns. The two cases of 2.5 and 5 //sec with their 
respective energy gains and pulse shapes are given because the first one equals the 
pulse length that was initially available on the SLAC accelerator and the second, in 
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conjunction with a 50 MW klystron, would satisfy the 50 GeV energy requirements of 
the SLC. The two repetition rates in combination with the corresponding pulse lengths 
require the same AC power. 

Note that while the SLED system is indeed a form of energy compression, it is 
not perfect because the cavities start out by rejecting the energy, then absorb some 
of it themselves, and finally do not emit nice rectangular pulses. Thus, for example, 
the power gain in the 5 psec case, instead of being 5 fzsec/0.8 /zsec, i.e., 6.25, is only 
3.15. The idea in itself, however, has been applied with considerable success and other 
laboratories such as CERN and the High Energy Institute at Beijing are adopting it 
for their linac designs. 

In the preceding pages, we have referred a number of times to the SLC or Linear 
Collider project at SLAC. Of all the known applications of electron and positron 
linacs, this is probably the most challenging because it puts to a test and stretches to a 
maximum the capabilities of these accelerators. For the SLC to become operational, we 
must improve our understanding and control of beam generation, injection, emittance, 
beam centering and focusing to keep the emittance from growing, energy spectrum, 
phase and amplitude stability of the klystrons, positron generation and many others. 
These problems are discussed elsewhere in this book by R. Stiening. 

The SLC will be the first test-bed for these ideas. Beyond this first 50 GeV-on-50 
GeV e ± linear collider, it is now being asked whether machines of this type could be 
built in the 1-2 TeV range. This question will of course be much easier to answer after 
some experience has been obtained with the SLC. However, one problem can already 
be focused on now. This is the problem of energy consumption. We will end these 
lectures with a simple example which might stimulate the reader to think further. 

It is not clear at this point where the ultimate limits of accelerating gradient will 
be nor what factors will determine the length of the accelerator. Suppose we want to 
build a 1 TeV-on-1 TeV e* collider. Say we have two choices of gradient: 20 MV/m 
(conservative, attainable now) or 100 MV/m (less conservative, yet to be tested). The 
respective lengths would be 50 km and 10 km/linac. Returning once more to Eq. (5.2) 
and rewriting it as 
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we see that even if we are able to bring the energy onto the axis of the accelerator 
at 100% efficiency, it will take a large number of joules to establish the fields. Table 
5 summarizes the results, assuming that r/Q = 5000fi/m at 2856 MHz. The power 
requirements are worked out for an example where the two linacs are pulsed 200 times 
per second to reach the overall desired luminosity. Each pulse is used to accelerate n 
bunches (4 < n < 10) which will at most remove 20% of the stored energy. Unless 
something very ingenious is done, the remaining 80% energy will be wasted. Note 
further that the real power consumption will be at least twice the rf power consumption 
because of all the inefficiencies involved in converting AC to rf. In a few years, at 
$0.10/kW-h, 1/2 GW of power will cost $50,000/hr. It is quickly seen that there will 
be a strong incentive to improve upon the numbers in this table. Two avenues suggest 
themselves. One is to increase the frequency since ur/Q scales as w 2 . The other is to 
invent a method of recovering the energy for some useful purpose. The challenge is 
yours! 

Table 5. 

rf Energy and Power Needs for a 1 TeV-on-1 TeV e ± Linear 
Collider Assuming No Losses in Establishing Fields 


Gradient (MV/m) 

20 

100 

Length/linac (km) 

50 

10 

rf Frequency (MHz) 

2856 

2856 

w s t (Joules/m) 

4.47 

111.73 

Total energy stored/linac/pulse (MJ) 

0.223 

1.117 

rf Power/2 linacs/200 pps (MW) 

89.36 

446.8 
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